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Abstract - A unified electro dynamic approach to the guided wave excitation by external 
sources in the waveguiding structures with bianisotropic media is developed. Effect of electric, 
magnetic, and magneto-electric losses in such media manifests itself in the special form of 
eigenmode orthogonality referred to as the quasi- orthogonality relation. It reflects the existence 
of the cross-power flow P^i and loss Qki for any pair (k, I) of modes which are rigidly linked to 
each other by this relation. The quasi-orthogonality relation remains true in the limiting case 
of lossless waveguides yielding the customary relations of orthogonality and normalization for 
propagating (active) modes and also their generalization for nonpropagating (reactive) modes. 

It is shown that the eigenmode set for a waveguiding structure is complete only outside 
the region of exciting sources. Inside this region the modal expansions of fields are incomplete 
and must be supplemented with the orthogonal complementary fields which extend the proper 
Hilbert space spanned by waveguide eigenfunctions. Among exciting sources there are the 
external bulk sources (currents, fields, and medium perturbations) and the external surface 
currents. Besides, the orthogonal complementary fields generate the effective surface currents 
on boundaries of the bulk exciting sources. 

The problem of waveguide excitation by external sources is solved by means of determining 
both the mode amplitudes for the modal field expansions and the orthogonal complementary 
fields inside the source region. The equations of mode excitation are derived on the basis of 
three approaches applying the direct use of Maxwell's equations, the electrodynamic analogy 
with the mathematical method of variation of constants, and the conjugate reciprocity theorem. 



CONTENTS 



1. INTRODUCTION 

2. GENERAL POWER-ENERGY RELATIONS OF ELECTRODYNAMICS FOR 
BIANISOTROPIC MEDIA 

2.1 Poynting's Theorem 

2.2 Mode Power Transmission and Dissipation 

3. ORTHOGONALITY AND QUASI-ORTHOGONALITY OF MODES IN LOSS- 
LESS AND LOSSY WAVEGUIDES 

3.1 Quasi-orthogonality Relation for Lossy Waveguides 

3.2 Mode Orthogonality in Lossless Waveguides 

3.2.1 Orthogonality and Normalization Relations for Active Modes 

3.2.2 Orthogonality and Normalization Relations for Reactive Modes 

3.3 Time-average Stored Energy for Active Modes in Lossless Waveguides 

4. ORTHOGONAL COMPLEMENTS AND EFFECTIVE SURFACE CURRENTS 
INSIDE SOURCE REGION 

4.1 Bulk and Surface Exciting Sources 

4.2 Orthogonal Complementary Fields and Effective Surface Currents 

5. EQUATIONS OF MODE EXCITATION 

5.1 Approach based on the Electrodynamic Method of Variation of Constants 

5.2 Approach Based on the Reciprocity Theorem 

5.2.1 Derivation of the Conjugate Reciprocity Theorem 

5.2.2 Derivation of the Equations of Mode Excitation 

6. CONCLUSION 
APPENDIX 
REFERENCES 



1 



1. INTRODUCTION 



Guided-wave electrodynamics in the modern understanding deals with the study of propaga- 
tion, radiation, excitation, and interaction of waves in a variety of waveguiding structures. The 
term "complex media waveguide" applied in the title of the paper implies the medium com- 
plexity of two types: 

(i) the physical complexity associated with medium properties diversified by the very nature 
(gas and solid state plasmas with drifting carriers; polarized and magnetized solids with differ- 
ent properties: piezoelectric, electrooptic, acoustooptic, magnetooptic, magnetoelastic; chiral, 
biisotropic, and bianisotropic media); 

(ii) the geometrical complexity due to using composite and multilayered structures. 
Electromagnetic theory has been developed up to its present state by extensive works and 

efforts of a great number of researchers and scientists. Besides pure scientific purposes, progress 
in classical electrodynamics at all stages of its advancement was encouraged by certain demands 
of technology. 

At the first stage such a stimulating factor was related to practical needs of then incipient 
radar and antenna engineering. The consequent experience on electrodynamic properties of 
mostly passive nondispersive media specified by phenomenological constants, which had been 
gathered over a number of years, was accumulated in many scientific publications. Among 
them we should refer, for instance, to such famous and popular books as [1-5] which now 
constitute the theoretical foundation of classical electrodynamics. Much attention was given 
to the study of electromagnetic properties of gas plasma as a medium for wave propagation. 
Later on the plasma wave aspects were extended to the behavior of charge carriers in solids 
considered as a solid-state plasma. At present the literature devoted to the electromagnetic 
properties of plasmas is immense and the following books [6-11] with their bibliographies can 
give a good indication of the scope of plasma electrodynamics. 

Another direction of electrodynamic aspects was inspired by developing the technology of 
microwave devices operating on wave principle. The first to be developed were vacuum devices 
using the space charge, cyclotron, and synchronous waves on an electron beam such as the 
traveling-wave tube, backward-wave tube, and others [12-14]. Later the similar idea to apply 
waves in solids for signal processing gave rise to new lines of solid-state electronics. They 
are due to applying the surface acoustic waves (SAW) in elastic piezo-dielectrics [15-17], the 
magnetostatic spin waves (MSW) in magnetized ferrites [18-23], and the space charge waves 



(SCW) in semiconductors with negative differential mobility of electrons [pX| , 24, |25|]. These 
waves refer to the quasistatic part of the electromagnetic spectrum of waveguiding structures for 
which a relevant potential field (electric for SAW and SCW or magnetic for MSW) dominates 
over its curl counterpart. This fact caused some electrodynamic formulations to be revised in 
order to separate such potential fields and take into account the space-dispersive properties of 



these media described by the proper equations of medium motion |15 ; £2], |2J, 25 1. 
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For the last decades the macroscopic electrodynamics of waveguiding structures has ex- 
perienced two powerful stimulating actions. The first is associated with needs of fiber and 
integrated optics and began about twenty five years ago. A number of theoretical propositions 
in electrodynamics were reformulated, as applied to optical waveguides, and have been em- 
bodied in devices. The literature devoted to this topic is enormous including the well-known 
books [26-32]. 

Nowadays we observe the renewed interest in electrodynamic problems caused by efforts to 
apply chiral, biisotropic, and bianisotropic media for the control of electromagnetic radiation in 
waveguiding structures. Phenomenon of optical activity in certain natural substances generated 
by their handedness property (chirality) was already known last century. The present revival 
of scientific and technological attention to this problem is inspired by the modern progress of 
material science and technology in synthesizing artificial composite media. Such media pos- 
sess unique properties to open new potential possibilities in their utilizing in optics and at 
microwaves. This has aroused a great wave of research followed by numerous publications, 
among them there are the general books [33-35] comprising bianisotropic issues and the spe- 
cial books [36-40] devoted entirely to this subject. 

Theoretical ground for many wave electrodynamics applications is the modal expansion 
method. In the case of the eigenmode excitation by external sources the question of com- 
pleteness of the eigenfunction basis chosen inside the source region is of crucial importance 
in practice. Unfortunately, most authors solve this question rather superficially assuming in- 
tuitively that the set of eigenfunctions found as the general solution to the boundary-value 
problem without sources is complete also inside the source region. However, this is not the case 
in general. 

From mathematical considerations given in Appendix A. 1 it follows that the above state- 
ment is valid only for the desired functions ip(x) tangential to the Hilbert space spanned by 
the eigenfunction basis {ipk(x)}- Generally, for most functions f(x) their series expansion in 
terms of the base functions (convergent in mean) is only a projection ip{x) of the function f(x) 
on the Hilbert space. In addition, there may exist a nonzero function c(x) orthogonal to this 
space, the so-called orthogonal complement, which in general must be added to the projection 
ip(x) in order for f(x) to be considered as the complete required function (see Eq. ( [A. 18 ) and 
relevant relations in Appendix A. 2). The above statement is fairly obvious for mathematicians 
but unfortunately was fully ignored in developing the modern topics of guided-wave electro- 
dynamics by most authors, not counting Vainshtein [0] and Felsen and Marcuvitz ||. Strange 
as it may seem, when developing the excitation theory of optical waveguides, many authors 
[26. 27, 29, [31], 32] have correctly applied the modal expansions to the transverse components of 
electromagnetic fields but entirely dropped the orthogonal complements due to the longitudinal 
exciting bulk currents. As will be shown, this causes the so-called effective surface currents 
to be lost. Similar situation also holds for electrodynamics of the waveguiding structures with 
chiral and bianisotropic media [34, 38 - 40] where the modal expansion method is practically 
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undeveloped and the problem of the orthogonal complements and effective surface sources, 
worked out below, is more complicated. 

The objective in writing this paper is to develop a unified electrodynamic theory of wave- 
guide excitation by external sources (bulk and surface) applicable equally for any media and 
waveguiding structures. Particular attention will be given to the study of the unexpandable 
orthogonal complements to the eigenmode expansions which should be expressed in terms of 
the given exciting currents as well as the desired mode amplitudes of the modal expansions. 
To this end, we begin with Sec. 2 devoted to deriving the basic energy-power relations of 
electrodynamics applied to the lossy bianisotropic media including Poynting's theorem in the 
differential and integral forms involving the self-power and cross-power quantities (flows and 
losses) transmitted and dissipated by the eigenmodes of a waveguiding structure. Sec. 3 deals 
with a generalization of the known orthogonality relation for the waveguides without losses to 
the so-called quasi- orthogonality relation for lossy waveguides which describes, as a special case, 
the orthogonality of the reactive (nonpropagating) modes in lossless waveguides. In addition, 
an expression for the time-average energy stored by the active (propagating) modes is proved. 
Sec. 4 is concerned with the consideration of external sources (currents, fields, and medium per- 
turbations) and electromagnetic fields inside the source region. The complete representation of 
the fields, besides their modal expansions, involves also the so-called orthogonal complementary 
fields which necessarily generate the effective surface currents. Sec. 5 contains two different ap- 
proaches to the derivation of the equations of mode excitation by external sources. The first 
approach applied only to the lossless waveguides is based on an electrodynamic analogy with 
the known mathematical statements such as the method of variation of constants and the rela- 
tions of functional analysis (see Appendix A). The second approach makes use of the reciprocity 
theorem in the complex-conjugate form to obtain the equations of mode excitation in the ge- 
neral form valid for both lossy and lossless waveguiding structures. Another alternative proof 
of the excitation equations for lossless waveguides starting directly from Maxwell's equations 
is adduced in Appendix B. 

In this paper we restrict our consideration to the case of time-dispersive media whose elec- 
trodynamic properties (isotropic, anisotropic, bianisotropic) are characterized by the frequency- 
dependent constitutive parameters considered as phenomenologically given. More complicated 
case of space-dispersive media such as elastic piezo-dielectrics, magnetized ferrites, nondege- 
nerate plasmas with drifting charge carriers whose electrodynamic description requires, besides 
Maxwell's equations, employing the proper equation of medium motion will be the subject of 
matter of the second part of the paper. 

In conclusion there are a few words concerning the notation applied: 

(i) tensors of rank (scalars), I (vectors), 2 (dyadics), and more than 2 (tensors) are de- 
noted as: A, A, A, and A, respectively; 

(ii) their products are denoted as: AB (for two scalars); A • B, A x B, and AB (for scalar, 
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vector, and dyad products of two vectors); AB • CD = AD(B • C), A • B = AijBjk, and 
A • B = AijkBkim (for scalar product of two vector dyads, dyadics, and tensors); AB : CD = 
(A • D)(B • C), A : B = AijBji, and A : B = AijkBkji (for double scalar product of two 
vector dyads, dyadics, and tensors). 

2. GENERAL POWER-ENERGY RELATIONS OF ELEC- 
TRODYNAMICS FOR BIANISOTROPIC MEDIA 

2.1 Poynting's Theorem 

In macroscopic electrodynamics, the electromagnetic properties of a medium are described 
by two field- intensity vectors, E (the electric field) and H (the magnetic field), and two flux- 
density vectors, D (the electric induction) and B (the magnetic induction), which are related 
by means of Maxwell's equations (written in the rationalized mks system): 

9B <9D 
VxE = - — , VxH = — + J, V-D = p, V-B = 0. (2.1) 

dt dt H y ' 

Mobile charge effects in the medium are specified by the charge and current densities p 
and J, whereas the bound charges arise as a result of polarization responses of the medium to 
electromagnetic actions characterized by the electric and magnetic polarization vectors P (the 
polarization vector) and M (the magnetization vector). These vectors yield the corresponding 
conrtibutions to the electric and magnetic inductions: 

D = e E + P and B = /z (H + M) . (2.2) 



The conventional procedure applied to Eqs. ( |2.1[ ) reduces to Poynting's theorem in the form 
involving the instantaneous values of power-energy quantities: 

_E + VS = -Ij - I P - I M (2.3) 

where w = (E • D + H ■ B)/2 is the electromagnetic energy density and S = E x H is the 
electromagnetic energy flux density (Poynting's vector). The terms on the right of Eq. ([2. 3D 



Ij=JE, (2.4) 



2 V dt dt 
lM = ^H-^r-^oM-— j (2.6) 

reflect specific properties of the medium under study and take into account the instantaneous 
power of interaction between the electromagnetic fields (E, H) and the charges - both mobile 
ones carrying the current J and bound ones generating the polarization P and magnetization M. 

In the literature the energy term dw/dt is conventionally identified with the sum E-<9D/<9t+ 
H-<9B/<9i, which is true only if D =e-E and B =/x-H where the tensors e and/x are symmetric 
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and time- independent. In this case only the first interaction term (2.4) is taken into account, 
whereas two others ( |2.5| ) and ( |2.6| ) are dropped without any justification. As will be evident 
from our subsequent examination including the second part of the paper, these terms play an 
important role in the power-energy theorem. 

For time-harmonic fields (with time dependence in the form of exp(iivt) ) one is usually 
interested in time-average values of the power-energy quantities denoted as (...). In this case 
(dw/dt) = so that Eq. ( p.3j ) takes the following form involving the time-average values of 
quantities: 

V-(S) = -(Ij)-(I p }-(l M ). (2.7) 

Below we concentrate on bianisotropic media for which there is no equation of motion. 
Their properties are usually described by the constitutive equations establishing macroscopic 
local relations among field vectors. It should be emphasized that magneto-electric effects 
(for instance, optical activity), by their microscopic nature, are brought about by nonlocality 
of polarization response on electromagnetic actions |53|, [jSj, [H], But their macroscopic 
manifestations are usually similar to those of actual time-dispersive media because for plane 
waves with the wave vector k = (u/c)n all the constitutive tensor parameters of such media 



become solely frequency-dependent (see Ref. plf ). 



There are a few forms of the constitutive relations for bianisotropic media [33 - 38] . Among 
them we choose the following form 

D = e E + | H, (2.8) 

B = C E + fi H, (2.9) 

as more convenient for our subsequent examination. 

Four constitutive medium parameters e, p,, £ and £ are considered as dyadic functions 
of frequency given phenomenologically. They comprise all special cases of the physical media 
without space dispersion: 

(i) for the isotropic medium 

e = el, p, = nl, | = C=0; (2.10) 

(ii) for the double anisotropic medium 

e^el, p^/il, | = C=0; (2.11) 

(iii) for the chiral (biisotropic) medium 

e = el, fi = fil, | = (x ~ *«) V e o A«o I, C = (x + V £ o Mo I, ( 2 -12) 

where x an d k are Tellegen's parameter of nonreciprocity and Pasteur's parameter of chirality, 
respectively [35, 38 1. It is known [33-35] that for a bianisotropic medium without losses the 
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dyadics e and /x are hermitian (self-adjoint) while the dyadics £ and £ are hermitian conjugate 
(mutually adjoint), that is 



■t 



/x = /x 



t 



(2.13) 



where superscript ' denotes transpose and complex conjugate (hermitian conjugate). Rela- 
tions ( p. 13 ) imply that in the general case of lossy media the antihermitian parts e a = (e— e ^)/2, 
p, a = (/x— /x t)/2 and the difference (£— C ) are responsible for losses (dielectric, magnetic, and 
magneto-electric, respectively). If the medium has also the electric losses related to its con- 
ductive properties and specified by the conductivity dyadic <x c , then in addition to Eqs. (|2[ 
and (l2.9|) there is another constitutive relation 



J = <7 r E . 



(2.14) 



Let us calculate the terms in the right-hand side of Eq. ( p. 7] ) by using their definitions (2.4) 
through (U) and the constitutive relations fl2j), p^9|) , and Q2.14|) : 



(Ij) = (J-E) = iRe{J.E*} = ^Re{E*-<r c -E} = : 



EE* 



1 /„ 0P „ 0E 
I P = ~(E - — -P- — 

V P/ 2 \ at dt 



(2.15) 
(2.16) 



I /e ~ - D ■ ?\ = - Re (e* • (iwD) - D* • (iwE)) 



at at / 4 

-Re{w^(V-e-E + E*-f -H-E-e*-E*-E-£*-H*)} 



-Re{icj[E*- (e-e f ) -E + 2E*-£-H 



| = iRe{io;(e a : EE*+|: HE*)}, 



V M/ 2 \ at ^ u at 



(2.17) 



= K H -1- B -^) = i^{h*. (w b)-b*.^h)} = 

= iRe{iw(H* -/x-H + H* -C-E-H-/i*-H* -H-C*-E*)} = 
= iRe{iw[H*- (/X-/X 1 ") •H-2E*-C t -H]} = ^Re{ico(fl a : HH* - ft : HE*) }. 
Substitution of Eqs. ( plBD - (gl7|) into Eq. Q 

gives the time-average Poynting theorem 

in the following form 

V- (S) + (q) = (2.18) 
involving the average Poynting vector 



(S) = - Re{E x H*} 



(2.19) 
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and the average power loss density 



(q) = <Ij) + (Ip) + (Im> 



= ^a e : EE* + ~<x m : HH* + ^Re{ef me : HE*} 
where we have introduced the total tensor of electric conductivity 



°"c + = °"c + iwe a = a c + iuj ■ 



e — e 



(2.20) 



(2.21) 



associated with conductor (cr c ) and dielectric (a^ = iue a ) losses of a medium, the tensor of 
magnetic conductivity 



i 



(2.22) 



associated with magnetic losses of a medium, and the tensor of magneto- electric conductivity 



a me = iu(£ - C 



?0J 



(f a + C a ) + (i h -c' 



(2.23) 



consisting of both antihermitian (with superscript a) and hermitian (with superscript h) parts 
of the cross susceptibilities £ and £. Unlike <x me , the dyadics <x e = <x c + iooe a and <r m = icj/i a 



are hermitian so that they produce the real (positive) definite quadratic forms in Eq. (2.20). 



2.2 Mode Power Transmission and Dissipation 

In order to obtain expressions for the power carried by modes along a waveguiding structure 
involving complex (anisotropic and bianisotropic) media and to find the dissipation of mode 



power it is necessary to go from the time-average Poynting theorem in differential form (|2.1§|) to 



its integral form. For this purpose let us integrate Eq. ( 2.18; ) over the composite (multilayered) 
cross section S = Yl^i formed from a few medium parts Si with interface contours Lj by using 
the two-dimensional divergence theorem (e.g., see Ref. §, p. 150) 





[ V-AdS = — ( z -AdS+ I n - A 

J Si OZ J Si J Li 



dl 



(2.24) 



where A is the arbitrary field vector and no is the outward unit vector normal to the contour 
Li and perpendicular to the longitudinal unit vector zq. 

Application of the integral relation ( 2.24} ) to V • ( S ) gives 

d 



V-(S)dS = ^ I z -{S)dS-J2f n+.(S+) + nr.<S- 

S OZ J S J Li 



di 



(2.25) 



where (S^) means values of the time-average Poynting vector taken at points of contour Li 
lying on its different sides marked by the inward (for either adjacent medium) unit vectors 
nf. The parts of interfaces between two adjacent nonconducting media do not contribute to 



the line integrals in Eq. (2.25) owing to continuity in tangential components of the electric and 
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magnetic fields. The only contribution may appear from the parts of Li due to conducting 
surfaces on which there is the known boundary condition [Q, ||, 35 1 



E T = Z s ■ (H T x m ) (2.26) 

where E T and H r are the electric and magnetic fields tangential to the surface and Z s is 
the surface impedance tensor. For the special case of the isotropic metallic surface with the 
conductivity a and the skin depth 5 = we have ||] 



Z s = (1 + i)U s I where TZ S = 1/aS = ^ivno/^a . (2.27) 

In this case the integrand of the line integral in Eq. ( p. 25 ) yields the surface loss power den- 
sity ( q' ) in addition to the bulk loss power density ( q ) entering into Poynting's theorem ( [2. 18 ). 

The result of integrating Eq. ( |2,18| ) over the cross section S of a waveguide and applying 
Eqs. ( 2.19[ ), ( |2.20| ), and ( |2.25 ) - ( |2.27| ) gives Poynting's theorem in the integral form 

dP / \ 

— + Q = 2.28 

dz 

where the total real power carried by electromagnetic fields in the direction of increasing coor- 
dinate z is equal to 

P = ( { S ) • z dS = - Re / (E x H*) • z dS (2.29) 
Js 2 Js 

and the total power loss per unit length caused by the bulk losses (obtained by integrating 

(q) over S = an d the surface (skin) losses (obtained by integrating (q') along 

L = Li) is equal to 

Q = + QW = [(q)dS + [ (q')dl = 

JS JL 

= \( (& e :BE*)dS+ \ f (<x m :HH*)dS + 
2 Js 2 Js 

+ i Re / (a me : HE*) dS + \ f K s (H r • H*.) dl . (2.30) 
2 Js 2 Jl 



Power relation ( 2.28 ) is valid only for the source- free region of a waveguiding structure 
whose electromagnetic fields can be expanded in terms of its eigenmodes (cf. Eqs. (|A.29|) ) 

E(r t ,z) = Y, A kMn)^ lkZ = 52a k (z)± k (r t ), (2.31) 

k k 

H(r t ,z) = J2 A k~Kk(r t )e-^ z = ^ a k (z) H fc (r t ) . (2.32) 

k k 

Every fcth mode is specified by the propagation constant j k = a k + if3 k and the eigen- 
functions {E^, H^} (where the hat sign over field vectors implies their dependence only on 
transverse coordinates r^, see Appendix A. 2), which are regarded as known quantities found 
from solving the appropriate boundary-value problem. The amplitudes A k are determined by 
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the exciting sources and called the excitation amplitudes. Inside the source region they depend 
on z as a result of source actions but for the source-free region A k (z) = const., as in the case of 
Eqs. ( [2.31 ) and ( |2.32| ). It is often convenient instead of A^z) to introduce the mode amplitudes 



a k (z) = A k (z)e-^ z 



(2.33) 



which take into account the total z-dependence related both to the mode propagation (exp(— 7^.2;)) 
and to the exciting sources (A k (z)), if any. 

Let us employ the modal expansions ( |2.31 ) and ( 2.32j ) to calculate the power flow P{z) and 



the power loss Q(z) given by Eqs. ( [2.29| ) and ( 2.30 ) for the source- free region. The final result 
of calculations is the following: 



P(z) = — / (E* x H + E x H*) • z dS = 

4 Js 



(2.34) 



Q(z) 



k I 



(E* • <r e • E) dS + 



k I 



(H* ■a m -Yl)dS + 



+ \ [ (E* - a me ■ H + H* • <4 e • E) dS + i / K s (H; • H T ) 

4 2 Jl 



dl 



1 



M ki aUz)a l (z) = J2Y,Qki(z 



k l k I 

where we have introduced the normalizing coefficients 



(2.35) 



N kl = (E^xHi+^x Hfe) • z dS (2.36) 
Js 

and the dissipative coefficients 

M kl = 2 / (E* k • <f e • E,) dS+ 2 / (H£ • a m • H,) dS + 

+ J s (K- °me • + H£ • <r ^ ■B l )dS+2j L TZ S (H* fc • H T , ; ) d2 (2.37) 

constructed of the cross-section eigenfield vectors (marked with the hat sign above them, see 
Eq. (pD). 

From Eqs. ( 2.36j ) and ( 2.37| ) it follows that the matrices {N k {\ and {Mjy} are hermitian, 
that is 

N kl = N? k and M kl = M* lk , (2.38) 

and have dimensions of watts and watts per meter, respectively, because the amplitudes A k 
and a k are dimensionless. 



The quantities P k (z) and Q k (z) appearing in Eqs. (|2.34|) and (|2.35|) for / = k in the following 
form 

Pk(z) = P kk (z) = ^N kk a* k (z)a k (z) = ^N k \A k \ 2 e" 2 ^ (2.39) 
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Qk(z) = Qkk(z) = \M kk al(z)a k {z) = i M k \A k \ 2 e~ 2akZ (2.40) 

are the real self powers transmitted and dissipated at point z by the kth mode which was 
excited at point z = with amplitude A k . 

Similarly, the quantities P k i(z) and Q k i{z) for I ^ k equal to 



P 



1,1 



(z) = -N kl al{z) ai {z) = lN kl AlA ie -^+^> (2.41) 



Q u (z) = \ M kl a* k (z) ai (z) = ~ M kl A\A X e" ^+^ z (2.42) 

can be interpreted as the complex cross powers transmitted and dissipated at point z jointly by 
the fcth and Zth modes which were excited at point z = with amplitudes A k and A\. Owing 
to ( p.38 ), the quantities defined by Eqs. ( |2.41| ) and ( |2.42j ) are also hermitian: 

P kl (z) = Pf k (z) and Qu{z) = Q* lk {z) • (2.43) 

From Eqs. ( |2.41| ) through ( p. 43 ) it follows that in a lossy waveguiding structure every pair 
of modes always transmits and dissipates the real combined cross powers 

Pfo(z) = P kl (z) + P lk {z) = 2ReP kl (z) = ^Re{N kl a* k {z)ai(z)}, (2.44) 

Q c kl {z) = Q k i(z) + Q lk (z) = 2ReQ kl (z) = i Re{M k i a k (z)ai{z)}. (2.45) 

Therefore, the double sums in Eqs. ( gJD and ( gTgg ) yield the real (time -average) total 
powers transmitted and dissipated by all modes in a lossy waveguide: 

p(z) = EE P ^) = E P ^) +EE p M = 

k I k k l^k 

= lj2 N kWk{z)? + iReEE^a^^^), (2.46) 

k k l>k 



c <z) 



Q( z ) = E E QM = E Qk(z) + EE<3« 

k I k k Ijik 

= -J^M k \a k {z)\ 2 + -ReEE^^)^), (2-47) 

k k l>k 

where N k = N kk and M k = M kk . 

In the next section we shall derive a relation named the quasi- orthogonality relation to link 
the cross powers P k i and Q k i for every pair of modes in a lossy waveguide or for every pair of 
the so-called twin- conjugate modes in a lossless waveguide. 
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3. ORTHOGONALITY AND QUASI-ORTHOGONALITY 
OF MODES IN LOSSLESS AND LOSSY WAVEGUIDES 



3.1 Quasi-orthogonality Relation for Lossy Waveguides 



Let us begin our examination with the general case of the composite (multilayered) wave- 
guiding structure containing bianisotropic media with bulk (electric, magnetic, magneto-electric) 
losses and surface (skin) losses. Consider the kth and Ith modes propagating in the source-free 
region of the waveguide which obey the curl Maxwell equations (|2.1| ) rewritten by using the 



(3.1) 
(3.2) 



constitutive relations fl2.8|) , (2.9), and ( 2.14 ) in the following form 



V x H fc(Z ) = (<r c + iuje ) • Efcjx) + iu>f • B. k ^ . 
A conventional procedure applied to Eqs. (|3.1| ) and (^) gives 

V.(Ejj;xH, + E,xH£) = 
= — 2 E£ • <x e • E; — 2 H£ • a m ■ Hi — ( E£ • <x me • + H£ • tx^e • E/) 



(3.3) 



where we have used formulas ( 2.21| ) - Q2.23 ). 

Application of the two-dimensional divergence theorem ( |2.24| ) to the left-hand side of 
Eq. ( |3.3| ), by analogy with formula ( 2.25| ) and by using the boundary condition ( |2.26| ), results 
in the following expression 





l fc • CT m • n; 



H,) dS 



— f (EJExH, + E,xHj).z d5 
= -2/ (Ej^.E,)^- 2 / (HK 

Representation of the fields for the kth and /th modes in the form 

E fc(0 (r tj z) = E fe(/) (r t )e-^) z , H fc(0 (r t , z) = H fc( ,)(r t ) e~W 



(3.4) 



(3.5) 



and their substitution into Eq. ( |3.4| ) give, by comparing with formulas fl2.36| ) and ( 2.37| ) for 
Nki and Mj~i, the desired relation 



(7fc+Tl)^fci= M kl 



(3.6) 



referred to as the quasi-orthogonality relation. It will play the same role in deriving the exci- 
tation equations (see Eq. ( |5.46| )) as the ordinary orthogonality relation. 

Expressions ( |2.41|) and ( 2.42Q relate the normalizing and dissipative coefficients Nki and 
Mki to the cross-power flow P^i and the cross-power loss Qki, respectively, transmitted and 
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dissipated jointly by the kth and Zth modes. The use of these expressions allows us to rewrite 
the quasi-orthogonality relation ([O]) in the power form 



(7* + li) Pki = Qki ■ (3.7) 

This formulation provides the following power interpretation of the mode quasi-ortho- 
gonality: outside the source region every pair of modes, independently of other modes, trans- 
mits the complex cross-power flow P k i rigidly coupled to the complex cross-power loss Q k i by 
the factor (7^ + 7;) consisting of the mode propagation constants ^yui) = a k(i) + *Afc(Z) > an d 
in doing so the combined cross powers P£> = P k i + Pik and = Q k \ + Qik always remain 
real. The quasi-orthogonality relation fl3.7|) means that outside the source region Poynting's 
theorem ( p. 28 ) takes place for any one of mode pairs (k, I): 

dPki 



dz 



+ Qki = 0. (3.8) 



Besides, every single mode has the real self-power flow P k = Pkk and self-power loss Q k = 
Q kk in the form of Eqs. ( gjg ) and Eq. flggt] ). These self powers are coupled to each other 



by the same relations ( |3.7[ ) and ( p.8| ) which for I = k yield the following expression for the 
attenuation constant: 

Qk M k 

ak= Jp k = ^Tk- (3 - 9) 

Hence, there occurs the following pattern of mode power transfer in the lossy waveguiding 
structures. Every fcth mode propagates from the source region with the fixed value of amplitude 
A k (the loss attenuation is taken into account by the amplitude constant a k appearing in 7^) 
which was excited by the sources. Outside them the mode, being a linearly independent 
solution to the boundary-value problem, does not interact with other modes owing to their 
linear independence. The kth. mode transfers the self power P k on its own and the cross powers 
Pki in conjunction with the other Zth modes which were also excited inside the source region 
and outside retain constant their excitation amplitudes Ai as well as the kth mode. 

3.2 Mode Orthogonality in Lossless Waveguides 

The orthogonality relation for a lossless waveguiding structure is obtained from the general 
relation ( |3.6| ) as the special case of M k \ = and has the following form 

(7fe+ ll)N kl =0. (3.10) 

In spite of the absence of dissipation, in the eigenmode spectrum of the lossless waveguide, 
besides propagating modes with a k = and 7^ = there are also modes having complex 
values of the propagation constant 7& = a/c + ifik with a k 7^ 0. These modes exist in the cutoff 
regime of propagation and their attenuation is of reactive (nondissipative) character associated 
with the storage of reactive power. As a token of this, it seems reasonable to refer to such 
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modes as the reactive modes to distinguish between them and the active (propagating) modes 
carrying an active (real) power. 

In the literature the reactive (in our terminology) modes are variously termed the complex, 
cutoff, nonpropagating, and evanescent modes. The last term is usually assigned only to cutoff 
modes with pure decay (a^ ^ 0) and without phase delay = 0). The latter feature of 
evanescent modes makes appropriate for them also the term "nonpropagating" because there 
is no phase propagation. But for the complex modes with otk 7^ and flk 7^ their reactive 
decay as exp(— a^z) is accompanied by the phase variation in accordance with the wave factor 
exp[i(ujt — Pk z )]- F° r this reason it is more preferable to refer to the complex modes as reactive 
modes rather than nonpropagating ones. However, we shall apply both terms, the reactive and 
nonpropagating modes, as well as their antitheses, the active and propagating modes, to reflect 
the fact that the former do not transfer any self power, whereas the latter carry it. 

Let us show that in any lossless waveguide, independently of its structure and media used, 
every reactive mode with number k has its own twin mode with number k (marked by tilde) 
so that their propagation constants are related by the equality 

T k = ~lk or a k = ~ a k , /% = Pk- (3.11) 



Such mode twins with pair of numbers (k, k) that satisfy the relation (3.11) will be referred 



to as the twin- conjugate modes. As is seen from Eq. ( |3.11 ), these modes have the same phase 



velocity = (5^) but decay in opposite directions (a^ = — ar). 

The existence of twin-conjugate modes possessing the property expressed by Eq. fl3.11| ) 
can be justified by means of the following reasoning. In our treatment of complex amplitude 
technique, we have chosen the wave factor in the form exp[z(u£ — k • r)] where k z = — 27. 
However, there is another alternative form exp[— i(cot — k • r)] with k z = — i^/ which differs from 
the first form in opposite sign of imaginary unity and having tilde above the wave vector. It is 
clearly evident that the alternative case can be obtained from our solution by applying complex 
conjugation, then 

k z = —ij = k* = ry* whence 7 = —7*. (3.12) 



Equalities for gammas in Eqs. ( |3.11 ) and ( p. 12 ) are fully coincident not counting different 



positions of tilde (the former marking the mode number k in subscripts will be used later 
on). This result substantiates the existence of twin-conjugate modes for which jk = — 7? 
or ^z,fc = k* r- In other words, any dispersion equation obtained as a result of solving the 
boundary- value problem for lossless systems has the complex roots with complex-conjugate 
values of the longitudinal wavenumber k z which appear in pairs. Such a pair of complex roots 
corresponds to the twin-conjugate modes. 

Sign of the amplitude constant can be used as the basis for classification of the reactive 
(nonpropagating) modes under two types, forward and backward, as is usually done for the ac- 
tive (propagating) modes but on the basis of a sign of the group velocity v grj k = [d(3k{u>) / du:]~ l . 
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In reference to the source region location between z = and z = L, all the modes (active and 
reactive) can be classified into two categories: 

(i) the forward modes marked by subscript k = + n > (active with v gr ^ n > or reactive 
with a +n > 0) which, being excited inside the source region, leave it (without or with reactive 
damping) across the right boundary and exist outside at z > L; 

(ii) the backward modes marked by subscript k = — n < (active with Vg rt — n < or reactive 
with a-n < 0) which, being excited inside the source region, leave it (without or with reactive 
damping) across the left boundary and exist outside at z < 0. 

3.2.1 Orthogonality and Normalization Relations for Active Modes 

The active (propagating) modes exist in the pass band of lossless waveguides where they 
have zero amplitude attenuation (a k m = 0), so that their propagation constants 7&(n = iflun 
are pure imaginary. In this case the orthogonality relation ( |3.10| ) rewritten in the form 

(J3 k - A) N kl = (3.13) 

along with expression ( [2.36 ) for N k i furnishes two alternatives: 



N H = [ (E£ x H|+E( x BE) • z dS = for l^k (3.14) 
JS 



or 



N k = N kk = 2Re { (Et x H fc ) • z dS / for I = k . (3.15) 
Js 

Expression ( |3.14j ) is the orthogonality relation between the different propagating modes for 
which j3 k — Pi ^ 0, whereas formula (3.15) defines the norm N k = N kk of the kth mode. It 



should be noted that Eq. ( |3.14| ) does not necessarily hold for different but degenerate modes 
with (3 k = Pi. In this case one can employ the conventional technique commonly used for usual 
waveguides 0, [| to ensure the orthogonality among degenerate modes by constructing from 
them such linear combinations that constitute a new orthogonal subset for which relation ( p. 14 ) 



is applicable. For this reason we shall no longer turn special attention to degenerate modes. 



From Eqs. ( 2.41 ) and ( p.l4[) it follows that two different propagating modes (with numbers 
k 7^ I) have zero cross-power flow (P k i = 0), i.e., they are orthogonal in power sense. Any 
mode carries along a waveguide only the self power P k defined by formula ( 2.3S| ), which gives 



the following power interpretation for the norm of active modes: Nk is equal to 4P£ where P£ 
means the time-average power carried in the positive z-direction by the kth mode with unit 
amplitude (\A k \ = 1). In some instances it may be more convenient to normalize the mode 



amplitude to unit power (\P k \ = 1 watt). Then \N k \ = 4 watts and according to Eq. ( (2.39Q 

P k = ±\a k \ 2 = ±\A k \ 2 (3.16) 

where subscripts should be read as k = ± n, with upper and lower signs corresponding to 
the forward and backward modes for which, respectively, N +n = 4 watts and N_ n = — 4 watts. 
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In the special case of a reciprocal waveguide wherein for every forward mode there is a backward 
one with the same law of dispersion, their norms are related to each other by the equality 



N +n = - iV_ n . (3.17) 

In conclusion, let us write the relation of orthonormalization for the active (propagating) 
modes in the following form 

r I for I zfi k , ,„ _ „. 

Nkl = NkkSkl = ]N kk = N k for l = k\ ( 3 - 1§ ) 



where the normalizing coefficient N k i and the norm N k are given by Eqs. ( 3.14] ) and ( p. 15 ), 
respectively. 

3.2.2 Orthogonality and Normalization Relations for Reactive Modes 

The reactive modes of a lossless waveguide are cutoff modes whose propagation constants 
lk(i) = a k(i) +iPk(l) are generally complex- valued or particularly real-valued for the evanescent 
modes with j3 k (i) = 0. So the general relation of orthogonality ( |3.10 ) holds for them and by 



using Eq. ( |3.11 ) for the twin-conjugate modes gives two alternatives: 

N H = [ (Et x H/ + E/ x Hi) • z dS = for I ^ k (3.19) 
Js 



or 



N k = N kk = f(E* k xHj + Ej xU* k )-z dS^0 for l = k. (3.20) 
J s 



Expression (3.1 9|) fulfils a role of the orthogonality relation for the reactive modes. As is 



seen from here, every reactive kth mode is orthogonal to all the Zth modes (reactive with ai ^ 
and active with a.\ = 0) for which 7/ + j k = 7/ — j k 7^ and iVjy = 0, including itself since 
7fe + 7fc = 2«fc 7^ and N kk = 0. The only mode nonorthogonal to the given kth mode is its 
own twin with number k for which 7^ + 7^ = 0. Formula ( 3.2Q ) defines the norm, N k = N kk 



for the reactive kth mode which is constructed of the fields of twin-conjudate modes (k, k). 
From Eqs. Q), fljgj] ), and (@) il follows that every reactive mode has no both the self 



power (P kk = 0) and the cross powers with the other modes (P k i = 0) for which 7^ 7^ — 7^, i. e., 
these modes are orthogonal in power sense. 

Each mode forming a twin-conjugate pair, being nonorthogonal to its twin, has its own 
norm defined by Eq. (|3~2Cp as 



N k = N kk or N~ k = N~ kk , (3.21) 
whence, according to the general property of hermitian symmetry for the normalizing coeffi- 



cients expressed by equality ( 2.38 ), it follows that 

N k = N* , (3.22) 
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i.e., the reactive twin-conjugate modes have the complex-conjugate norms. 

Although the reactive mode has no self power (Pkk = 0), the twin-conjugate modes in pair 
carry the real combined cross power (cf. Eq. I l| 

P k = P kk + P kk = 2 P k~k = 1** {N k A* k A~ k } ee \ Re {N~ k A\A k ] (3.23) 



where subscripts should be read as k = + n and k = — n. This is a consequence of relations ( p. 11 ) 



for the twin-conjugate modes and the definition of forward and backward reactive modes: if 
the feth mode is a forward one with N k = N +n = N+ n - n , then the kth mode is a backward 
one with N- k = iV_ n = AU,+n = N$ n . 

As evident from Eq. ( |3.23| ), to transfer the real power by reactive modes it is necessary that 
both constituents of a twin-conjugate pair should have nonzero amplitudes (A k and A k ) and to 
be in such a phase relationship that their combined cross power P?r would be other than zero. 
Similar situation usually takes place in the regular waveguide of finite length bounded by two 
irregularities and excited at frequences below its cutoff frequency [1-3]. Reflections from these 
irregularities can form inside this length two evanescent (cutoff) modes with numbers k = +n 
(forward mode) and k = —n (backward mode) constituting the twin-conjugate pair for which 
[3± n = and 7 +n = a +n = — a_ n = — 7* n . It is easy to see that the norms for the forward and 
backward evanescent modes are pure imaginary- valued and related to each other by the general 
relation ( |3.22| ). Superposition of fields for the two evanescent modes with opposite decay sense 
furnishes nonzero real cross-power flow along a short length of the cutoff waveguide. 

In conclusion, let us write the relation of orthonormalization for the reactive (nonpropagat- 
ing) modes in the following form 



where the normalizing coefficient N k i and the norm N k are given by Eqs. ( gls| ) and ( pp|) 



respectively. From comparison of Eqs. (13111 ) and flpg ) it is seen that the latter relation is of 



general form because it comprises the former one for the active modes as a special case obtained 
by replacing subscript k with k so that, in particular, the norm N k = N kk takes the form given 
by Eq. fTJ). 

It is pertinent to note that all the above expressions for the norms and the relations of 
orthogonality and orthonormalization can contain the total field vectors in place of their cross 
section parts related to each other by Eqs. (|3l]|), i.e., the hat sign over the field vectors can 



be dropped. This is obvious for the active modes and follows from the equality 7fc + 7| = for 
the reactive twin-conjugate modes. 

If in a waveguiding structure there are both the active (propagating) and reactive (nonprop- 



agating) modes, the total power flow (|2.46[) carried by them, in accordance with the aforesaid, 
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is given by the following expression 

P = Pact + Preact = ^2 Pk + ^2 ^kk 

k k 
(active) (reactive) 

= \Y, N k \a k (z)\ 2 + ±ReJ2' N k a* k (z)a~ k (z) (3.25) 

k k 
(active) (reactive) 

where prime on the sum sign means summation of the twin-conjugate modes in pairs rather 
than that of the single reactive modes. 

3.3 Time-average Stored Energy for Active Modes in Lossless Waveguides 

The time-average Poynting theorem written in the form of Eq. ( 2.1§| ) for time-harmonic 
fields does not contain a stored energy density. In order to find it one usually applies variational 
technique (e. g., see Ref. (^9|). To this end, it is necessary to obtain a relation between variations 
of the electromagnetic fields (5E k , 5H k ) for the A;th mode and perturbations of the frequency 
and medium parameters (5(uje), 5(u)p,), 5(tu£), 6(uC)) which bring about these variations. 



The kth mode is governed by Maxwell's equations ( |3.1| ) and (3.2) with a c = for a lossless 
medium whose other parameters satisfy the requirements ( p. 13 ). By taking variations in these 
equations we obtain 

V x 5E k = - iuip, ■ 5H k - iuC, ■ 5E k - i5(ujfi) • H fc - iS(uC) ■ E k , (3.26) 
Vx<5H fc = iuoe ■ 5E k + iu£ ■ 5H k + i5{ue) ■ E k + i8(u£) • H k . (3.27) 
A conventional procedure applied to Eqs. ( |3.26| ) and ( 3.27 ) reduces to the following relation 



V • (E* k x 5H k + 5E k x H*] 



iuj 



Ej-(e-et).<JE fc . 



+ HI ■ [fi ■ 5U k + Efc • (| -C f ) • SU k + U* k -(C -V) ■ SE k 



-i E* k - S(ue) ■ E k + H* k ■ S(ufi) • H fc + E^ • 6{u>£) ■ H fc + H* k • 6(u£) ■ E k 



where the terms inside the first square brackets vanish because of relations ( |2.13| ) for lossless 
media so that 

V • (E* k x 5H k + 6E k x H* k ) = 



8(ue) : E k E* k + 8(ujfi) : H k H* k + 2Re{c«) : H fc E fc 



(3.28) 



Electromagnetic fields of a propagating mode and their variations can be written on the 
basis of Eq. fl3,5|) as 



E* e- i(3kZ 



H k = H k e- l ^ z 



5E k = {5E k -t5p k z-E k )e- i ^ z , 
5H k = (5H k -i8(3 k zIl k )e- i ^ z . 



(3.29) 
(3.30) 
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When substituting Eqs. ( 3.29 ) and ( |3.30| ) into Eq. ( 3.2S ) and applying the integral rela- 



tion ( |2.25| ) where the line integral vanishes owing to continuity in tangential components of the 
fields, the integration over the cross section S yields 

Sp k ^Ref (B k xHl)-z dS = 
2 Js 

= - [ \ 5{ue ) : E fe E£ + 8{ufi) : H fc H£ + 2 Re{%f) : H fc E£} 1 dS. (3.31) 

The left-hand side of Eq. fl3.31|) involves the time-average power flow P k as a multiplier of 
5(3 k - By using the known relation 

Pk = v gr ,k W k where v gr , k = f ~^j~ ) ( 3 - 32 ) 

is the group velocity, we obtain from Eq. ( |3.31| ) the desired expression for the time-average 
energy stored per unit length of a waveguide (dropping the mode index k): 

w = \ /(e*.%^.e + h*.^2.h) dS + 

+ ^Re / (e* ■^^-■U)dS= [ wdS. (3.33) 



2 Js\ duj / js 



In accordance with expression ( |3.33 ), the time-average stored energy density for a lossless 



bianisotropic medium is given by the following formula 



1 



4 



d(uje)/duj d(uo^)/duo 



™ = 7 E* H* • • „ (3.34) 



d(LoC,)/doo d(ujii)/duj 




which is a generalization of the usual expression for lossless anisotropic media J7|, [KJ, 22]. 



4. ORTHOGONAL COMPLEMENTS AND EFFECTIVE 
SURFACE CURRENTS INSIDE SOURCE REGION 

4.1 Bulk and Surface Exciting Sources 

Up to the present, the external sources exciting the composite (multilayered) waveguiding 
structures involving isotropic, anisotropic, and bianisotropic media have been dropped. From 
this point onward, the special attention will be given to investigating the behavior of modes 
inside the source region. In doing so, we assume that all the eigenfields {E^H^,} in the form 
of Eq. ([^), including their eigenfunctions of cross-section coordinates {E/%,Hfc} (marked by 
hat over them) and their eigenvalues of propagation constants jk = a k + iPk, are known from 
solving the corresponding boundary- value problem. As shown in Appendix A, these eigenfields 
constitute an infinite countable set of the vector functions quadratically integrable on the cross 
section S of a waveguiding structure. This set can be taken as a basis of the proper Hilbert 
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space to expand the required fields E and H not only outside sources, as was done by Eqs. ( [2,31 ) 



and ( p.32 ), but also inside the region of external sources. In general, this eigenvector basis is 



not complete inside the source region since it cannot take into account entirely the potential 
fields of the sources. This requires to supplement the modal expansions E a and H a with un- 
known modal amplitudes Af-(z) by the orthogonal complements Ej, and Hf, (see Eqs. ( [A. 27 ) 



and (A.2Sp). Hence, the desired issues to be obtained inside the source region are both the 



longitudinal dependence of modal amplitudes and the orthogonal complements to the modal 
expansions. 

In the most general case there exist three physical reasons to excite the waveguiding struc- 
ture under examination: 

(a) the external currents - electric 3 e ext and magnetic J™. 4 , 

(b) the external fields - electric E ea; f and magnetic H e xt > 

(c) the external perturbations of bianisotropic medium parameters Ae, Ap,, A£, and A£. 
Owing to these medium perturbations, the total electric (E+E ea:f ) and magnetic (H+H ex t) 

fields create the excess electric AD and magnetic AB inductions linked by the constitutive 
relations Q2.8| ) and (fT9|), that is 

AD = Ae • (E + B ext ) + A| • (H + H ext ) , (4.1) 
AB = AC • (E + B ext ) + Ap ■ (H + U ext ) . (4.2) 

These excess inductions bring about the induced displacement currents — electric J| n( ^ — 
iujAT) and magnetic J™ d = iuA'B which, being added to the external currents 3% xt and J^, 
yield the bulk exciting currents 

3% = r ext + 3\ nd = 3% xt + iuAY) , 3? = J£t + JZa = J£t + iwAB (4.3) 
entering into the curl Maxwell equations ( |2.1| ) in the following form 

V x E = - iioB - 3^, (4.4) 

V x H = icuD + 3 e b . (4.5) 



The conduction current J of a conductive medium defined by Eq. (2.14) is now assumed to 
be incorporated with the electric displacement current iuD, whereas the induction vectors D 
and B are taken, as before, to be related to the intrinsic electromagnetic fields E and H inside 



the medium in question by the same constitutive relations ( |2.8| ) and (2.9). So the permittivity 



tensor e is now regarded as a sum (e + a c /ioj) whose antihermitian part defines the total tensor 
of electric conducttivity <r e =<r c + given by Eq. ( 2.21| ). 



Besides the bulk exciting currents and 3V 1 , there may exist the surface exciting currents 
3 e s and J™ which give discontinuities of the appropriate tangential components of fields at 
points of the surface whereon these sources are located, written in the form of the following 
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boundary conditions: 

n+ x E+ + n s 7 x E" = - J™, (4.6) 

n+ x H+ + n s x H = J e s . (4.7) 

Here the field vectors with superscripts ^ mean their values taken at points of the source 
location contour L s lying on its different sides marked by the inward (for either adjacent 
medium) unit vectors nf. 

4.2 Orthogonal Complementary Fields and Effective Surface Currents 

The general electrodynamic eigenmode treatment (see Appendix A. 2) based on the well- 
known mathematical formulations (see Appendix A. 1) yields the complete representation of 
the desired field vector F inside sources as a sum of the the modal expansion giving a 
projection of F onto the Hilbert space and the complement C orthogonal to the Hilbert space 
sec Eq. flOp ). Thus, the elect romagnetic fields inside the source region have the complete 
representation given by Eqs. flA~27| ) and flO§), namely (cf. Eqs. JOlD and (|2~32j)): 

E(r t , z) = E a (r t , z) + E 6 (r t , z) = ^ a k{z) E fc (r t ) + E 6 (r tj z) , (4.8) 

k 

H(r t ,z) =U a (r t ,z) +H b (r t ,z) = £ a k {z) K k (r t ) + U b (r t , z) , (4.9) 

k 

where E b (rt,z) and H b (rt,z) are the required orthogonal complements. The mode amplitude 
a k( z ) = Ak(z) exp(— 7^2;) allows for the total dependence on z due to both the unperturbed 
propagation of the A;th mode with a constant 7^ and the perturbed amplitude A k (z) as a result 
of source actions. Using Eqs. ( |2.8| ) and Q2.9| ) gives the similar expressions for the induction 
vectors 

D(r t ,z)=D (r t ,*) + D 6 (r t ,;s) = £ a k {z) D fc (r t ) + D 6 (r t , z), (4.10) 

k 

B(r t , 2) = B (r t , «) + B 6 (r t , z) = £ a fc (z) B fc (r t ) + B 6 (r t , z). (4.11) 

In order to find the orthogonal complements let us substitute Eqs. ( [4.8] ) through ( |4.11| ) into 
Maxwell's equations ( |4.4j ) and ( [4. 5| ) with taking into account the fact that the eigenfields (|3.5|) 
satisfy the homogeneous (with no sources) Maxwell equations ( |3.1| ) and ( |3.2| ). Some transfor- 
mations yield 

dAi 

^^(z xE fc ) = -VxE fc - iuB b - , (4.12) 

k 



dA 



^(z xH,) = -VxH t + iwD fe + 3 e b . (4.13) 

The left-hand side of Eqs. ( [4.12 ) and (4.13) has only transverse components. From here it 
necessarily follows that there must exist nonzero orthogonal complementary fields. Otherwise 
(when Ef, = H b = = B& = 0) these equations become physically contradictory because 
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then they require the longitudinal components of the arbitrary bulk currents Jf and J5 71 to be 
always equal to zero, which is of course not the case. 

Hence, the required orthogonal complementary fields should be chosen so as to make the 
longitudinal component of the right-hand part of Eqs. ( |4.12 ) and (4.13) vanish, that is 



V • (z x E 6 ) - z • (i^Bft + Jft" 

V • (z x H b ) + z • (iu~D b + 3 e b 



0. 

0, 



(4.14) 
(4.15) 



where the identity zo • (V x a) = — V • (zo x a) has been used. 

Since the field parts Eft and Hft form the orthogonal complement to the Hilbert space 
spanned by the base eigenvectors {Efc,Hfc}, they must be orthogonal to the fields of any 
eigenmode in power sense given by the relation similar to Eqs. ( |3,14j ) and ( |3,19| ) (cf. Eq. ( |A.26[ )): 



/ (E£ x Hft + Eft x H£) ■ z dS = 
Js 



(z x Eft) • U* k - (z x Hft) ■ 



dS = . 



(4.16) 



In virtue of arbitrary choice of the kth eigenmode taken from the base set, zero equality in 



Eq. ( 4.16| ) can occur if and only if 

z x Eft = 
z x Hft = 



or 



or 



Eft = zo Ef, , 
Hft = zq -ffft , 



(4.17) 
(4.18) 



i.e., both complementary fields are longitudinal. In order for their magnitude to be found, it 
is necessary to insert Eqs. (|2~8j) , (|2~9|) , (|4.17|) , and (|4.18j ) into Eqs. (|4.14|) and fl4.15|) , then 

1 



z o • Dft = e zz Eh + £ zz Ht, 



J bz > 



1 



z o • Bft = Czz Eb + fi zz H b — — — Jbl ■ 

IUJ 



(4.19) 
(4.20) 



From here it finally follows that 



E b 



J bz v J bz 

iooe zz (1 — v e v m ) 

J bz v J bz 
iiJfJ,zz (1 - v e v m ) 



where we have denoted 



( zz C : z z 



izz i ■ zozo 



(4.21) 
(4.22) 

(4.23) 



t-zz e : zozo fJ-zz M : z z 

Hence, both the orthogonal complementary fields are longitudinal and produced by the 
longitudinal components of the bulk exciting currents. 
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Existence of the complementary fields E b and immediately reduces to appearance of the 
so-called effective surface currents i e s e j and J™ e j- Consider the bulk source region having the 
cross section Sb with a boundary contour L b and write the complete electric field inside and 
outside this area: 



E(r t ,z) 



J2 k A k (z)E k (r t ,z) + E b (r t ,z) 
J2 k A k (z)E k (r t ,z) 



E - inside S& , 
E + - outside S b . 



(4.24) 



Analogous expressions can be written for the magnetic field H(r^, z). 

The eigenfields E k and H k , being obtained for the situation without sources, are generally 
continuous at points of the line Lb- Then the tangential components of the complete fields E 
and H prove discontinuous: 



n^" x E + + n b x E = — rif, x E5 , 
x H + + x H = - n b x H 6 , 



(4.25) 
(4.26) 



where = ri£ = — is the outward unit vector normal to both the line L b and the longi 
tudinal unit vector zo- The comparison of these relations with the boundary conditions 
and (4.7) yields the desired effective surface currents 



's,ef 



n b x H 6 



J ™e/ = n b x E fc 



iuJUzz (1 - v e v m ) 
JUL b ) - v m J%{L h ) 



iuje 7 



(4.27) 



(4.28) 



where J bz {L b ) and JJ^(L b ) mean the longitudinal components of the bulk currents taken at 
points lying on the boundary L b of their existence area S b and r = zq x n^is the unit vector 
tangential to the contour L b . 



The general expressions ( [4.21 ) and ( 4.22 ) for the complementary fields E& and H b and the 
general expressions ( 4.27|) and ( 4.28j) for the effective surface currents J e s e j and J™ e / take the 
following simplified form in special cases of: 

(i) the isotropic medium with parameters ( 2.1C| ) {y e = u m = 0) 



E,, 



zp 

iuje 



J, 



bz 



fin 

J s,ef 



J bz( L b) 



and 
and 



H,, 



J s,ef 



J, 



T 



bz > 



(ii) the anisotropic medium with parameters ( 2.11 ) (u e 

zo 



E,, 



J, 



's,ef 



iu>e z 



Jbz{L b 



and 
and 



H 6 
J 



v m = 0) 

Zf) 



Jbz(L b 



jrn 



bz ■> 



s,ef 



Jbz(L b 



(4.29) 
(4.30) 

(4.31) 
(4.32) 
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Formulas (pL29D and (^3l|) are in agreement with those obtained first by Vainshtein and 



Felsen and Marcuvitz Q, respectively. As for the effective surface currents ( 4.30| ) and ( [4. 32 ) 



Vainshtein did not consider them at all but the excitation integrals in the theory of Felsen and 
Marcuvitz allow for them implicitly, which will be shown later (see Sec. 5.1). 

Therefore, in the absence of medium bianisotropy the bulk currents, electric J| and magnetic 
J™, generate the effective surface currents, respectively, magnetic J™ e j and electric J^e/- As 
is seen from Eqs. ( 4.27] ) and ( [4.28 ), the bianisotropic properties of a medium intermix the 



contributions from the bulk currents into the effective surface currents owing to the longitudinal 
components £ zz and C, zz . 

The newly obtained effective surface currents 3 e s e j and J™ e j, as well as the actual surface 



currents 3 e s and J™ entering into the boundary conditions (4.6) and (4/7), make contributions 
to the excitation amplitudes A k along with the bulk currents and J™. 

The next step should be done toward deriving the differential equations to find the functions 
A k {z) inside the region of bulk and surface sources. For this purpose we shall apply three 
independent approaches set forth in the next section and Appendix B. 

5. EQUATIONS OF MODE EXCITATION 

5.1 Approach based on the Electrodynamic Method of Variation of Con- 
stants 

Mathematical method of variation of constants is applied to solve an inhomogeneous dif- 
ferential equation (with driving terms) by representing its general solution in the form of a 
superposition of the known linearly independent solutions of the proper homogeneous equation 
with coefficients which are no longer considered constant and assumed to be the desired func- 
tions of an independent variable J42|. Electrodynamic analog of the mathematical method of 
variation of constants is built by representing the fields E(r^, z) and H(r^, z) inside the source 
region in the form of the expansions (^^) and ( |4.9[) in terms of eigenfunctions of the proper 
homogeneous boundary- value problem (without sources) whose amplitude coefficients A k {z) 
are the desired functions of z rather than constants, as they are outside sources. 

According to the conventional mathematical technique, the method of variation of constants 
is to give differential equations for the mode amplitudes in the following form 

^- = fk(z), k = l,2,... (5.1) 

where the functions fk(z) take into account the longitudinal distribution of exciting sources 
(bulk and surface). Integration of Eq. (|5.1|) yields the required dependence 



A k (z) = A° k + J f k (z)dz = A% + AA k (z) . (5.2) 

The integration constant A^ should be determined from a boundary condition given at one 
of two boundaries (z = or z = L) of the source region depending on the type of modes for a 
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lossless waveguiding structure: 

(i) for the forward modes (active and reactive, k = + n) supplied at the left input 



A +n (0) ± or A +n (0) = , (5.3) 
(ii) for the backward modes (active and reactive, k = — n) supplied at the right input 

A-n(L) / or A- n (L) = . (5.4) 



Substitution of Eq. ( |5,2p into Eqs. (O) and ([Oj) allows us to represent the complete 
solution for the electromagnetic fields inside the source region as the sum of the general solution 
(E gen , Hgen) to the homogeneous boundary- value problem (without exciting sources) involving 



the constant amplitude coefficients and the particular solution (E 
inhomogeneous problem (with exciting sources), namely: 

E(r t , z) = E 9en + E par = 

= Y, A °k B k(r t , z) + Yj A Mz) Efe(r t , z) + B b (r t , z) 



H r 



par) *-*-par 



to the proper 



H(r t ,z) 



par 



A° k H fc (r t , z) + A M?) H fc (r t , z) + H 6 (r <s z) 



Thus, the general technique of solving the electrodynamic problem of waveguide excitation 
by external sources based on the method of variation of constants gives rise to the representation 
of the desired electromagnetic fields as the sum of the general and particular solutions adopted 
in the theory of linear differential equations. The next task is to obtain a specific form for the 
excitation equation like Eq. ( |5.1| ). 

To derive the equation of mode excitation let us vector- multiply both sides of Eqs. (JO]) and 
(|4.9|) by and — Ef, respectively, and add them. Then after scalar-multiplying the result 
of summation by zo and integrating over the cross section S of a waveguide we obtain 

J (E* x H + E x H*) -z dS = J (E/* x H fc + E fc x H z *) • z dS + 



(5.5) 



The last integral in the right-hand side of Eq. (|5.5| ) vanishes because of the orthogonality 
relation ( |4.16| ) or ( |A.26| ). According to the orthonormalization relations (|3.18| ) and ( |3.24 ), the 
integral under the sign of summation is equal to NiSik for the active (propagating) modes and to 
NiSr, for the reactive (nonpropagating) modes. Hence, from Eq. (5^) we obtain (cf. Eq. ( A.. 31 )): 



+ / (E,* x U b + E h x H*) • z dS 



(i) for the active modes (with replacing subscripts / — > k) 

1 



E x Hfc) • z dS = A k e 



-if3 k z 
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or 



A k = ^- f (E* k x H + E x Hj£) • z dS 



(5.6) 



where the norm N k is defined by formula ( |3.15| ) , 

(ii) for the reactive modes (with replacing subscripts I — > k and I — ► fc) 



or 



— y s (E| xH + ExH!)-z d5 = A fc e" 
A k = ±- J s (El xU + B x Hi) ■ z dS 



(5.7) 



where subscript k corresponds to the mode which together with the kth mode constitute the 
twin- conjugate pair and have the propagation constant 7? = — 7^ and the norm N? = N k 
defined by formula ( |3. 20 ). From comparison of Eqs. ( |5.6| ) and ( |5.7| ) it is seen that the latter 
expression can be considered as the general form valid not only for the reactive modes but also 
for the active modes with replacing k by k. 

It is pertinent to note that expression ( [5.6D for a k is in agreement with Eq. QA..31 ) obtained 
by minimizing the mean-square difference D n (defined by Eq. ( |A.33|) ) between the mode se- 
ries expansion \1/ and the partial sum S72 of the nth. order (given by Eq. ( |A.32|) ) to provide 
convergence in mean for the modal expansion. 



Formulas (5.6) and (5.7) give a rule to find the mode excitation amplitude if the electro- 
magnetic fields are known. However, this is usually not the case because the exciting currents 
(bulk and surface) are assumed to be given rather than the fields. In order to go from the fields 
to the currents, let us differentiate the general relation ( |5.7| ) with respect to z: 

dA k _ d 

dz 



dz 



|(E|xH + ExH|)-z c/S 



(E! xH 



E x Hi) dS + 



+ E 



n, 



(E~ x H + E x H!) + +n:.(E|xH + ExH 



dl 



(5.8) 



where the last equality is written by using the relation similar to Eq. (2.25). 

The complete fields E and H inside the source region satisfy the inhomogeneous Maxwell 
equations ( [4.4| ) and ( [4.5D , whereas the fields E~ and H| of the kth mode obey the following 
homogeneous equations 



VxE! 



V x Hi 



iojB t 



k 



By using the constitutive relations (gj), ( J2T9[ ) and Eqs. (gj), (g75|), (|5 
easy to prove that 



(5.9) 
(5.10) 

and ( 5,10| ) it is 



V • (E| x H + E x H|) 



(J^-Ei+J^-Hi) 



(5.11) 
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(e -et) : EE| + (ft : HH| + (| f ) : HE? + (C -ft) : EH~ _ 

where the square bracket equals zero for a lossless medium owing to Eq. ( gig ). 

The contour integrals in the right-hand side of Eq. (fT^) include two contributions: 

(i) from the actual surface currents 3 e s and J™ which are located on a contour L s and meet 
the boundary conditions ( |4.6|) and Q4.7|), 

(ii) from the effective surface currents i e s e j and J™ e y given by Eqs. ( p~27D and ( p~28l) which 
are located on a contour bounding the bulk current area S& and meet the boundary condi- 
tions ( p5| ) and ( p6|) . 

On the strength of the aforesaid we can write 



E 



+ 



n+ • (E| x H + E x H!)+ + nr . (E| x H + E x H|) 

n+ • (Ef x H + E x HS)+ + nj • (E| x H + E x H|)' 

n+ • (E| x H + E x H|)+ + n 6 ~ • (E| x H + E x H?)- 

= - / (JI • E| + J™ • H|) dl - (J e s>ef ■ E| + J™ e/ • H| 

Eqs. ( |5.8| ), ( 5.11 ), and ( |5.12| ) finally give the desired equations written as 
(i) for the excitation amplitudes A^^z), k = ±n: 

dA k 



dl = 
dl + 
dl = 
dl. 



(5.12) 



dz 



i-/ (Jf.E| + Jr-H|)d5- 



Nr 



(^•E| + J--H!)d/ 



(J 



(5.13) 



(ii) for the mode amplitudes Ofc(z) = Ak{z) exp(— jk z )i k = ±n: 



dak 

—r- + Ik ak 

dz 



±- j a (Jg-Ei- J II - 



k JS b 



Nr 



1 



— /_ (j; • E| + J- • HI) dl- — /_ (J« e/ • E| + J- , • H|) efl 



(5.14) 



The excitation equations ( |5.13j ) and (5.14) written for the amplitudes of reactive modes 
hold true also for an active mode if one assumes k = k and 7^ = i(3k- The excitation integrals 
in the right-hand side of these equations represent the complex power of interaction between 
the external currents (bulk and surface) and the eigenfields of the /cth mode (for active ones) 
or those of its twin-conjugate kth. mode (for reactive ones). 

As distinct from the theory developed, Vainshtein fully excluded from consideration the 
reactive (nonpropagating) modes and the effective surface currents and restricted his analysis 
only to the reciprocal waveguides with isotropic media. In this case every forward-propagating 
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mode (k = + n) has a backward counterpart {k = — n) of the same type so that their common 
norm is defined by Vainshtein as 



N n = (E +n x H_ n - E_ n x H +n ) • z dS . 
Js 



Unlike the definition ( 3.15| ), Vainshtein's norm has no power sense and does not allow a 
generalization to nonreciprocal waveguides to be made. 

Theory of Felsen and Marcuvitz [||, unlike Vainshtein's theory, takes into consideration 
anisotropic (not bianisotropic) media but also does not allow for the reactive modes. The 
excitation integral in their equation similar to our Eq. ( 5. 14] ) has a visually different form which 
does not involve the effective surface currents explicitly. In order for their implicit existence to 
be displayed, let us convert our excitation integral containing the bulk currents. 

To this end, it is necessary to transform the products of longitudinal components such as 
J£ z E kz and H kz (where k = k for reactives modes) . The use of the constitutive relations 
D = e • E and B = fi ■ H for a double- anisotropic medium in Eqs. Q5,9j ) and ( |5.10| ) written for 
the kth mode gives the longitudinal projections of these equations: 

(V x E* k ) z = (V t x E* kt ) ■ z = uj (v* zx H* kx + lf zy H* ky + ^ zz H* kz ) , (5.15) 

(V x Hfe) 2 = (V t x H* kt ) • z = - iu (e* zx E* kx + e* zy E* ky + e* zz E kz ) . (5.16) 

Taking into account that for a lossless medium e*j = €ji and = on the basis of 
Eqs. (5.15) and ( |5,16j ) we can obtain the following expressions 

je F * _ 1 m v tt* \ T e e ^ E kx + e VZ E ky f 
J bz ^kz - - -■ \ V t X ™-kt) ' J bz J l 



iuje 7 

1 



and 



bz — 

^•^tl-^'tJ' 1 *-^ (5 - 17> 

jm tt* 1 (T7 v p.* \ T m ^ H kx + HyzHjy jm _ 
J bz M kz - ~ l V * X &kt) " J bz J bz - 

lUJflzz Hzz 

= L Vt .f KtX ^k\ + L( Vtx m.\ . E * _ th*j r (5.i8) 

ILO \ fi zz / ILU\ [l zz ) fX zz 

where following [|| we have introduced the auxiliary vectors 

e tz = x e xz + y e yz and fi tz = x fi xz + y ii yz ■ (5-19) 

It is easy to see that the terms in Eqs. ( |5.17| ) and ( 5.18| ) containing the transverse divergence 
operator ■ , after integrating over the bulk current area Sf,, yield the following results 

- / V t • (Kit x^-)dS = -I n b - U* kt x^)dl^ 

IU JS b \ fJ-zzJ IU JL b V HzzJ 

JlefKdl, (5.20) 



28 



s b 



l kt 



bz 1 dS 



1 



n 6 - H* w x-^ 

Li, \ c 2z 



dl 



's,ef 



n* k di, 



(5.21) 



where expressions ( 4.32 ) for the effective surface currents have been used. 

As is quite evident, the terms ( 5.20[ ) and ( |5.21| ), being inserted in the excitation integral with 
the bulk currents by means of equalities ( 5.17| ) and ( |5.18| ), fully compensate for the contribution 
from the excitation integral with the effective surface currents entering into Eq. ( 5.13j ) and 
(5.14). Then the excitation equation ( ^.14 ) written for propagating modes takes the form 
entirely coincident with that of Felsen and Marcuvitz H (in different notation): 



—j h 7fc «fc 

dz 



dS- 



where following M we have introduced the effective bulk currents 



1 



l^zz 



J| t +-(V t x-^l - —J§ z , 



€ l z 



(5.22) 



(5.23) 



J b,ef 



Jm 
bt 



\Hz 

l^zz 



J bz > 



(5.24) 



with the transverse vectors e^ z and fi tz being defined by formulas ( 5.19 ). The currents ( |5.23 ) 



and ( 5.24 ) were introduced by Felsen and Marcuvitz in different designations but of the same 
structure. 

being formed 



From Eq. ( |5.22D it follows that the effective bulk currents ^% e f and J™ e y, 



as mixtures of the longitudinal and transverse components of the actual electric and magnetic 
currents, interact only with the transverse eigenfield components E^ t and H^ t of the kih mode, 
but in doing so take into account the contribution from the effective surface currents i e s e j and 
J^ e j defined by Eqs. fl4.32| ). It should be mentioned that the contour integral in Eq. fl5.22| ) 
allowing for the contribution from the actual surface currents i e s and J™ is absent in the 
appropriate equation of Felsen and Marcuvitz |8| . 

In view of fundamental importance of the mode excitation equations, Appendix B displays 
another derivation for the lossless waveguding structures based on the direct use of Eqs. (|1|) 
and ( 4.13j ) which are an exact consequence of Maxwell's equations ( |4.4[) and (|4.5| ). The general 
case of lossy waveguides is studied below on the basis of the reciprocity theorem in complex- 
conjugate form. 
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5.2 Approach Based on the Reciprocity Theorem 

5.2.1 Derivation of the Conjugate Reciprocity Theorem 

The basis of deriving the reciprocity theorem in complex-conjugate form is constituted by 
two systems of Maxwell's equations like Eqs. (|4.4|) and Q4.5]): 



VxEi = - iwBi - , 



V x Hi = iuDi + JSi , 



VxE;= iu;B* 2 - J 



62 > 



V x H; = - iujT>% + J 



62' 



(5.25) 
(5.26) 



written for two different electromagnetic processes (marked with subscripts 1 and 2) excited 
by different external currents (bulk and surface), with the frequency and the constitutive pa- 
rameters of a waveguiding medium entering into relations fl2.8| ) and ( p.9| ) assumed to be the 
same. 



Application of the conventional technique to Eqs. ( 5.25| ) and ( |5.26 ) yields 



V- (Ei xHj + EjX Hi) = — (J e bl ■ E* + Jl* ■ EO - (J£ • H* + J^* • Hi) 
- iJ(Di-E;-D^-Ei) + (Bi H^-B^ Hi) 



(5.27) 



After transformation with using the constitutive relations (|2.8| ) and (|2.9| ) the last term in 



ILO 



the right-hand side of Eq. ( 5.27 ) accepts the following form 

'(Di-ES-DS-Ex) + (Bi-H^-B^-Hi) 
(e-et) :EiE^ + (fi-fi^-.Kjl^ + (£-C f ) : + : 



In accordance with the aforesaid in Sec. 4.1, the permittivity tensor e is regarded here 
as a sum (e + a c /iuj) so that its antihermitian part determines the total tensor of electric 
conductivity <x e = <x c + cr^ taking into account both dielectric {(td) and conductor (<r c ) losses 
of a medium. Magnetic losses (<r m ) are taken into account by the antihermitian part of the 
permeability tensor fi, whereas the tensor <x me = iui 

(£ - C f ) reflects the magneto-electric 
losses due to bianisotropic properties of a medium. The use of Eqs. ( ggjp - ( gjgjp converts 
relation (5.27) into the differential form of the conjugate reciprocity theorem 



V • S12 + q\2 



,(b) 
12 



where we have denoted 



Si 2 =e 1 xh; + e;xh 1 , 

912 = 2 (<r e : E1E2 + <r m : H1H2) + (<T me : H1E2 + it^ : E1H2) , 



,(&) 
12 



l J M ' ^2 + J 62 - -klJ _ l J M ' M 2 + J b2 ' 



(5.28) 

(5.29) 
(5.30) 

(5.31) 
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Superscript (b) reflects belonging an appropriate quantity to bulk properties of a system, while 
the surface properties will be marked by superscript (s). 

To obtain the integral form of the reciprocity theorem it is necessary to integrate Eq. fl5.28|) 
over the cross section S of a waveguiding structure with using the integral relation similar to 
Eq. ( 2.251 ) which involves the contour integrals taking into account two physical phenomena: 

(i) the skin losses expressed by the boundary condition ( |2.26| ) with the surface impedance 
tensor ( |2.27l ) given along a contour L , 

(ii) the discontinuity in tangential components of the fields caused both by the actual surface 



currents i e s and J™ located on a contour L s with the boundary conditions fl4.6|) and fl4.7|) and 
by the effective surface currents J e s e j and J™ e f located on a contour L& with the boundary 
conditions ( |4.25| ) and ( [4.261 ). 

For the sake of brevity it is convenient to write both surface currents as the overall surface 
sources 



located along the combined contour L s = L s + L&. 

Substitution of Eq. (|5.29|) into the integral relation ( 2.25| ) yields 

d r / „ \ 



(5.32) 





Is 


V • S 








-/ 


n i 


h -(Ex 


JL+L S 






+ f n b - 


(Eal 







dl + 



dl 



^J s (Vi x H*+E* x Hx) -z dS + 2jT (H T i-H* 2 )dl + 
r sl • E* + J™ • H*) + (J- • Ej + JS* • Hi) 



+ 



dZ + 



+ 



Lb 



(Js,e/1 " E a 2 + J™e/1 " H. 



u2/ 



+ ( J s*e/2 ' E al + J ™e/2 ' H al) 



dl. 



Here we have used: (a) the boundary condition ( |2.26| ) on the contour L with surface 
impedance ( |2.27l ), (b) the boundary conditions (fO|) and ( [4.7D with the actual surface currents 
•'si (2) anc ^ *'si(2) gi ven on the contour L s , (c) the effective surface currents J e s e ji( 2 ) = — n 6 x 
Hji( 2 ) and J™ e /i(2) = n b x ^bl(2) defined on the contour Lj,. Therefore, the line integrals in 
the previous formula yield two resulting contributions: 
(i) from the skin losses on the contour L of a conducting surface 



Q12 — 2H S (H r i • H* 2 ) , 
(ii) from the overall surface currents on the contour L E = L s + 



(5.33) 



» 

12 



( J El - ^2 + J E2 ' -kl) _ ( J El - M 2 + J E2 ' M l 



(5.34) 
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Therefore, the reciprocity theorem in the integral form is given by the relation 

dP 12 (z) 



dz 



+ Qia(z) = Ru(z) 



(5.35) 



where we have introduced the following integral quantities (complex-valued) (cf. Eqs. ( [2.29D 
and flpO|)) 



Pn{z)= I S l2 (r t , z) ■ z dS = [ (Ei x + x Hi) • z dS , 
Js Js 

Qn(z) = Q§(*) + qW(z) = f qi2 ( rt ,z)dS + [ q' 12 (r t ,z)dl = 

JS JL 

= 2 / (<x e : E1E2) dS + 2 / (<r m : H X H^ dS + 

Js Js 

+ [ (* me : H1E2 + <fL e : EiH5) dS + 2 / K s (H rl • H* 2 ) , 
Js Jl 

R 12 {z) = R? 2 \z)+r[ s 2 \z) = [ r^{v u z)dS + / r${v u z)dl 

Js h JLr. 



(5.36) 



(5.37) 



Sb J Lx, 

(J e bl • E2 + • K* 2 ) + (J 6 1 • Ei + J%? ■ Hi) 

l J El ' ^2 + J El - ±1 2J + 1 J e2 • ^1 + J e2 - H l/ 



(5.38) 



It is easy to see that with no sources (when R\ 2 = 0) the second system (with subscript 2) 
of Maxwell's equations ( |5.25| ) and ( |5.26| ) describes the same fields as the first (marked by 
subscript 1) only with taking complex conjugation. This makes it possible to replace subscript 2 
with 1 and what is more to drop them. In this case the integral reciprocity theorem ( |5.35| ) 
turns into the integral Poynting theorem ( |2.2&j ) in which the real power flow P and the real 
power loss (bulk and surface) Q per unit length of a waveguide are equal to 



P = 7 / Sh • z dS = ( ( S ) • z dS , 
4 Js Js 

Q = \ J s mdS + IJjudi = J s (q)dS + J L (q')di, 



(5.39) 
(5.40) 



where their expressions in terms of fields are given by Eq. ( |2~29l) and ( p0|) 



5.2.2 Derivation of the Equations of Mode Excitation 

Inside the source region the reciprocity theorem in the integral form (|5.35| ) is the basis for 
obtaining the excitation equations. To this end, the fields marked by subscript 1 (which will 
be dropped for the exciting currents) are assumed to be the desired fields excited by the bulk 
and surface sources (J b [ = J^' m 7^ and = J E ' m / 0) and represented in the form of 
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expressions ( [4.8D and ( |4.9| ) (with replacing summation index k by I), whereas those marked by 

re,n 

ha 



subscript 2 are the known fields of the kth. mode outside the source region (J^™ = J^™ = 



given in the form of Eq. (|3.5| ) . 

Substitution of Eqs. (gj) and (g§ into Eqs. ( ggg ), ( ^37|) , and flOgj ) yields the following 
expressions: 

Pi fc (^) = / Si fc (r t5 z) • z = / (Ei x H* k + E* k x • z dS = 
Js Js 

= J2N k iMz)e-^ + ^ z , (5.41) 



Qik(z) = / qik(r t ,z)dS + / q' lk (r t ,z)dl = 
Js Jl 

i 

R lk (z) = flg(z) + Hg>(z) ^ / r g(r t ,z)dS + / r«(r tl 



(5.42) 



z)dl 



Rf\z)e-% z + R k s '(z)e-% z = R k {z)e 



?(')/ 



7k* 



(5.43) 



where the normalizing and dissipative coefficients N k i and M^; have the form of Eqs. ( |2.36| ) and 
(|2 . 37p and the quantity Rk{z) = R k (z)-\-R k (z) consists of two exciting integrals (bulk and sur- 
face) : 

-/ (jt--E% + J?.H£)dS, (5.44) 

** Oh 



Rt\z) 



Rl s) (z) 



J s ' E fc + • Hfc ) dl 



's,ef ' ^k 



Et + 3T P f • Ht I dS . 



(5.45) 



These integrals involve the cross-section eigenfield vectors (marked with hat) and their 
dependence on z is due to that of the external currents 3 e b ,m {z) and J^ m (z). 



It should be pointed out that Eqs. (|5.4l[) and ( p.42[) come only from the field contributions 
of the mode expansions E a and H a since the orthogonal complementary fields E& and H5, 
being proportional to the longitudinal component of external currents, do not contribute into 
P\ k and cannot influence the intrinsic losses in a medium related to Q\ k . 

Substitution of Eqs. ( |5.41 ), ( 5.42 ), and ( 5.43| ) into the integral reciprocity theorem ( |5.35 ) 
(with replacing 2 by k) gives a relation 



J2\ N u^r ~ [dt + li)N kl -M kl 



Me 



Rk = R^ + R k ^- 



(5.46) 



The quasi-orthogonality relation of the general form ( |3.6| ) make the square bracket in 
Eq. fl5.46|) vanish so that it reduces to the desired set of the excitation equations written 
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in the following form: 

(i) for the excitation amplitudes A\(z) 



5>H^^e"^= Rf{z) +R { k s \z), k = l,2,... (5.47) 



(ii) for the mode amplitudes ai(z) = Ai(z)e 7(2 

= Rt\z) + R ( k s \z), A; =1,2,... (5.48) 



i 



dai(z) , . 



dz 



Discussion of the excitation equations obtained will be put off until the similar equations 
for the waveguiding structures with space-dispersive media are developed in the second part of 
the paper. 

Up to this point the waveguiding structures under study are assumed to be closed with a 
screening metallic boundary, whose spectrum of eigenmodes is always discrete. 

In conclusion, it is pertinent to show features of the excitation theory peculiar to open 
waveguiding structures (without losses) in which an outside homogeneous medium extends to 
infinity in one or both transverse directions. As is known p7| , 29, |3l[], for the open waveguides 
in addition to the discrete part of the spectrum of bound modes (with the outside medium 
fields localized near outer boundaries of the waveguiding layer), there is a continuous part of 
the spectrum related to radiation modes (with the fields extending to infinity in the outside 
medium). Unlike the eigenfields Efc(rj , z) and Hfc(r t , z) of discrete modes marked by the integer- 
valued subscript fc = l, 2, . . . and expressed by Eqs. (|3.5|), the fields of a radiation mode 



E(r t ,z;k?) =E(r t ;k?)e-W) 2 , H (r t , z; k?) = H(r t ; k?) e" i/3(k ° )z (5.49) 

are specified by the transverse wave vector k° = xo k° + yo k° of the outside passive medium. 

In this case the modal expansions of the fields E a (rt, z) and H a (rt, z) inside the source re- 
gion, besides the series expansion in terms of discrete modes, involve also the integral expansion 
in terms of radiation modes (cf. Eqs. (|4.Sj ) and (|4.G|)): 



E«(r t ,*) = 5> fc (z)E,(r t )e-^ + f A(z; k°) E(r t ; k?) e" l/3 ^ )z dk° , (5.50) 
k J 

H a (r t ,z) = J2M^k(r t )e-^ z + f A(z; k?) H(r t ; k?) e" l/3(k ° )z dk° , (5.51) 

k J 

where integrating over k° and k° is taken along the real axes from — oo to oo. 

The orthonormalization relation for radiation modes can be written by analogy with rela- 
tion ( |3.18 ) for discrete modes in the following form 

iV(k°, kf)= / [ (E*(r t ; k?) x H(r t ; k?') + E(r ( ; k?') x H*(r t ; k?) | • z (i5 = 

= 6(k t -k° t ')N(k°) (5.52) 
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where the Dirac delta function <5(k° — k°') = 5(k° — k°')5(k° — replaces the Kronecker 
delta function S^. It should be mentioned that since 5(k° — k°') has dimensions of (length) 2 , 
the dimensionality of the norm N(k°) and the excitation amplitude A{z; k°) for the radiation 
modes is equal to watts/m 2 and m 2 , respectively, as distinct from the bounded modes for 
which the similar quantities are taken in watts and as dimensionless. 

The equation for the excitation amplitude A(z; k°) of the radiation mode has the form 
similar to Eq. ( 5.13j ): 



dA(k?) 



/ (Jg-E*(k?) + J^H*(k°))dS- 



dz NQtfiJs^ 

- ^yj^(j£-E*(k?) + J?.lT(k?))cfl (5.53) 

where the coordinate variables are dropped for simplicity. 

The orthogonal complementary fields E& = zqE^ and = zqH^ obtained in the form of 
Eqs. ( 4.21 ) and fl4.22j ) remain valid for open waveguides. 



6. CONCLUSION 

We have shown a unified treatment of the electro dynamic theory of the guided wave excitation 
by external sources applied to any waveguiding structure involving the complex media with 
bianisotropic properties. Allowing for losses in such media has reduced to the power loss density 
in Poynting's theorem due to the magneto- electric conductivity a me defined by Eq. fl2.23j ), in 
addition to the usual electric and magnetic conductivities. 

Application of the desired field expansions in terms of eigenmode fields gives the self-power 
and cross-power quantities (flows and losses) transmitted and dissipated by the eigenmodes 
of a lossy waveguide, as well as the time-average energy density stored by the propagating 
modes in a lossless waveguide which involves the additional contributions from bianisotropic 
properties of a medium. 

The basis of developing the excitation theory for lossy waveguides is the novel relation ( |3.6| ) 
called the quasi- orthogonality relation whose general form is always true including the propa- 
gating (active) and nonpropagating (reactive) modes in lossless waveguides considered as a 
special case. Among the external sources exciting the waveguiding structure we have included 
the bulk sources (currents, fields, and medium perturbations) and the actual surface currents. 
Inside the source region the modal expansions ( |2 .31 ) and ( j2.32|) have proved to be incomplete 



and must be supplemented with the orthogonal complementary fields ( |4.21 ) and fl4.22| ), as it 



is done by Eqs. ( [Q| ) and (fL9| ). Generally these complementary fields generate the effective 
surface currents (14,27] ) and fl4.28|) . So in the most general case the external source region 
contains the bulk currents J^'"\ the actual surface currents Jf' m , and the effective surface 
currents J e s '™f brought about by the longitudinal components of the bulk currents. 
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The equations of mode excitation in the form of ( 5.47| ) or ( 5.48| ) have been derived by using 
three approaches based on: (i) the direct derivation from Maxwell's equations (see Appendix B), 
(ii) the electrodynamic analogy with the mathematical method of variation of constants (see 
Sec. 5.1), (hi) the reciprocity theorem in the complex-conjugate form (see Sec. 5.2). 

All the results obtained are valid for the time-dispersive media specified by macroscopically- 
local and frequency-dependent parameters. An extension of the theory to space-dispersive me- 
dia which require for their description the special equations of motion with regard for nonlocal 
effects will be examined in the second part of the paper where the orthogonal complementary 
fields are explained as a part of the contribution from the potential fields of external sources. 



APPENDIX 

A. BASIC RELATIONS OF FUNCTIONAL ANALYSIS AND 
THEIR ELECTRODYNAMIC ANALOGS 

A.l Mathematical formulation (in notation of [42]) 



Unlike [42], we shall examine the general case of nonorthogonal base functions which gives 
the orthogonal basis special case. 

Consider a countable set of complex functions ipi(x), tp2(x), ■ ■ ■ quadratically integrable in 
the sense of Lebesgue on a given set S of points (x). The class L2(S) of such functions (regarded 
as vectors) constitutes an infinite-dimensional unitary functional (vector) space if, in addition 
to two binary operations of the vector sum ipk(x) + ipi(x) and the product akipk(x) by a complex 
scalar a&, one defines the inner product of ipk{x) an d ipi{x) as 

(ipk,ipl) = \ ^* k {x)^{x)^i{x) dx (A.l) 
Js 

where the weighting function 7(x) is a given real nonnegative function quadratically integrable 
on S, in particular, may be j(x) = 1. 

If Gram's determinant det[(ipk, ipi)] built up on the inner products of the form (A.l) differs 



from zero, the functions ipk( x ), k = 1, 2, . . . are linearly independent in L2 and can be chosen as 
a basis of the unitary functional space, with their mutual orthogonality not being necessarily 
required in general. The given set of functions ipk(x) spans a linear manifold comprising all 
linear combinations of ip\ (x) , ip2 (x) , . . . . 

Let us compose a partial sum of the nth order 

n 

s n (x) = Y / a k n) Mx) (A.2) 

k=l 

with scalar coefficients ajf* not yet defined. Given a function ip(x) fully belonging to the linear 
manifold spanned by ip\ (x) , ip2 (x) , • • • , these coefficients can be found from the requirement 
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that the weighted mean-square difference 



D r 



■y(x)\s n (x) - il)(x)\ 2 dx 



(A.3) 



between s n (x) and ip(x) would be minimum. With the help of Eq. ( A.2 ) the quantity D n can 
be rewritten in the following form 



Dn= /[f>l n) >£(*)-^*(*) 
JS 1 k=l 



7(x) 



^a^ipiix) - V(i 
1=1 



dx . 



(n) 

Then the conditions of its minimality with respect to the set of coefficients a k are written as 



(A.4) 



dD ri 



da 



(n)* 



J2a[ n) ipi{x) -ip(o 
l=i 



dx = 



d 2 D r , 



da^da™ Js 



if;l(x)j(x)^k(x) dx > . 



(A.5) 



(A.6) 



Eq. ( |A.q ) complies with the requirement of quadratic integrability initially imposed on the 
base functions ipk(x), while the condition ( A.5 ) yields the following system of equations to 
find aj^ : 



Y^Nmu^ =R k , k = l, 2,... 
i=i 



where we have denoted 



(A.7) 

(A.8) 
(A.9) 



Metric convergence in L2 is defined as convergence in mean (with index 2) of the sequence of 
partial sums s n (x) (with coefficients a^ from Eqs. ( |A.7| )) to the function ip(x), i. e., s n (x) 
ip{x) as n^oo, which occurs if and only if 



N ki = / iPl( x h(x)^i(x) dx = (ipk,ipl) , 
Js 

Rk= ipt(x)j(x)ip(x) dx = (ipk,ip) ■ 



j(x)\s n (x) — ip(x)\ 2 dx —* as n — > 00 



(A.10) 



Using Eqs. (|A.7|) through (|A.9|) and the equality N M = N* k allows Eq. (|A.4[) to take the 
following form 

n n n n . 

D n = E E ^4 n N w - E ^4 n) * - E ^ (n) + / 7(*)i^)i 2 dx = 

k=ll=l k=l 1=1 JS 

n n „ 

= -EE^W+ / i(x)mx)\ 2 dx. 
k=n=i s 

From here for limiting case ( |A.1Q ), when a k n \z) —xik(z) as n— >oo, it follows that 

„ 00 00 

(i/>,M = / r(xh(xMx)dx = Y / Y, N u a *k a i • ( A - n ) 



k=l 1=1 
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Relation (A. 11) is realizable only for functions ip(x) quadratically integrable on S (with 
the weighting function j(x)), i.e., for which there exists an integral on the left. This relation 
expresses completeness of the set of the base functions tpk{ x ) (also quadratically integrable 
on S) inside the class of functions ip(x). The completeness property establishes the space L2 
as the Hilbert space for which the series expansion 



00 



^) m = n £«^) (A.12) 
fc=l 

interpreted in the sense of convergence in mean given by formula ( [A. 10 ) is valid. Uniqueness 
of this expansion arises from the following reasoning. 

By conradiction, let two different series expansions J2k ^k^kix) and Sfc a fcV ; fe( a; ) corre- 
spond to the same function ip{x) in the sense of convergence in mean. To determine the 
expansion coefficients a' k and a k ' there are two systems of form ( A. 7 ) with the same right-hand 
sides Rk- When resulted from them, the difference system of equations Y^i^kiWi — a") = 
0, k = 1, 2, . . . gives a[ = a'{ by virtue of iVjy / 0, i. e., the initial series expansions coincide. 
If on the contrary one assumes that the same series expansion J2 k akipk(x) corresponds to two 
different functions ip'(x) and ip"(x), then the difference function ip~(x) = ip'(x) — ip"{x) has 
the expansion coefficients identically equal to zero. So the right-hand side of the complete- 
ness relation ( |A.ll ) vanishes, which necessarily provides ip~(x) = 0, i.e., the initial functions 
coincide. 

The completeness relation ( A.ll| ) is a generalization of the conventional Parseval identity 
(see Eq. ( [OBI) ) to the case of nonorthogonal bases. All the aforestated convince us that linear 
independence and completeness of the set of base functions i/jk(x) are fundamental properties 
of the basis, whereas their mutual orthogonality is not obligatory requirement and merely 
facilitates the problem of finding the expansion coefficients ak- Indeed, for the orthogonal basis 



so that 



ty k ,ik)= il>ttx)7(x)i/> l (x)dx = for k^l (A.13) 
Js 

Nki = N k 5 k i with N k = f ^U x h{x)Mx)dx=\\^k\\ 2 (A.14) 



where \\tpk\\ = \/ C0fcj ipk) = V~Nk is conventionally called the norm of a function tpk( x ) |Hl- I n 
addition, we extend this term to quantities Nki recognizing the self norm Nk = Nkk for I = k 
and the cross norm Nki for I 7^ k. 



Hence, in the special case of the orthogonal basis satisfying Eq. ( A.14 ): 
(i) the system of coupled equations (|A.7|) falls apart into separate equations yielding 



a k = ^ = = i- / ^(x) 7 (*Ms) dx , (A.15) 



N k (^k,tpk) N k Js 

(ii) the general completeness relation ( A.ll| ) gives the conventional Parseval identity 

„ 00 
{^A)= / r{xh(x)ij(x)dx = Y,Nk\ak\ 2 • (A.16) 
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It should be remembered that the use of the known Gram-Schmidt orthogonalization pro- 
cess [42 1, in principle, allows one to construct the orthonormal basis. 



The above completeness property of a basis expressed by relation ( |A.ll ) or ( A.16| ) concerns 



only such functions ip{x) that fully belong to the linear manifold spanned by the functions 
ipi(x), tp2( x ), ■ ■ ■ ■ However, for the most general functions /(x) this is not the case. 

Any given function f(x) quadratically integrable on S (with the weighting function j(x), 
in general) can formally be associated with the function ip(x) represented by series ( |A.12| ) if 
one assumes that its coefficients at satisfying Eqs. ( |A.7| ) through ([A.9|) are due to f(x) and 
not to ip(x), i.e., the quantities Rk(x) contain f(x) in place of tp(x) under the integral sign 
of Eq. ( [A.g| ), Let us prove that the difference c(x) = f(x) — ip{x) is orthogonal to every base 
function ipk(%) ha the sense of relation ( A.13| ): 



(ipk,c) = (ipk, f~tp) = (ipk, f) ~ (i>k,ip) = 

= W>fc, f) ~ ^{ipk,in)ai =Rk~Yl N ki a i = ( A - 17 ) 
I i 

where the relations iVjy = (ipk,ipi) and Rk = (ipk,f) have been used. 

Thus, any arbitrary function f(x) not belonging fully to the Hilbert space (spanned, for 

instance, by eighenfunctions of a boundary- value problem) can be represented in the following 

form 

f{x) = V(x) +c(x) = J2 a k Mx) + c (x) . (A.18) 

k 

Here the function ip(x) written as a series expansion in terms of base functions (convergent 
in mean) and considered as tangential to the given Hilbert space is called the projection of fix) 
on this space, while c (x) is a function orthogonal to the Hilbert space and reffered to as the 
orthogonal complement because (ipk, c) = 0. For such a function f(x) instead of the generalized 



Parseval identity (A. 11) there exists the generalized Bessel inequality 

„ oo oo oo 

(/,/) = ||V + c|| 2 = / f*(xh(x)fix)dx>J2J2^kia* k a l or >^iV fc |a fc | 2 (A.19) 

S k=l 1=1 k=l 

where the last single sum corresponds to the orthogonal basis. 



A. 2 Electrodynamic treatment (in notation of [8]) 

Let us consider a relevant aspect of the electrodynamic modal theory on the basis of analogy 
with the foregoing mathematical relations. 

Given an infinite set of eigenfunctions of a boundary-value problem defined on the cross 
section S of a waveguiding structure with the transverse radius vector rt and the longitudinal 
axis z, any eigenfunction ^(r^) and its adjoint (hermitian conjugate) ^t(r^) are denoted in 
the two-vector notation as 

**(r*) = fl*^] and *t(r t ) = (Ej(r t ) H* fe (r t )) (A.20) 
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where the hat over field vectors means the absence of their dependence on z. 

By analogy with Eq. ( |A.1| ), the inner product of two eigenfunctions *&k( r t) and ^(r^) can 
be defined in the following form 



(* fc> *,) = / *J(r t ) • f • *,(r t ) dS (A.21) 
J s 



with the weighting function given in the form of a special dyadic 



r=( °T "T 1 ) (A-22) 



z x I 

where zq is the unit vector of the axis z and I is the unit dyadic such that (zo x I) • a = 
— a • (I x zo) = zo x a for any vector a. The weighting dyadic T is constructed so as to make 

the double scalar product *| • f ■ * z = ■ ff ■ *jj = • f ) • *, = f : under the 
integral sign of Eq. ( |A.21[ ) be equal to (E£ x Hj + E; x H£) • z . 

Thus, the cross norm Nki for the A:th and Ith. modes and the self norm = N^k for the 
fcth mode, according to Eqs. (A. 8) and ( |A.21| ), can be represented as 

iV fc , = (* fc ,* { )= / *l-r-*id5 = 

J s 

= [ (El x Hj + E; x H£) • z (A.23) 
Js 

and 

iV fc = (* fe ,* fe )= / - f -* fe d5 = 

= 2 Re / (E£ x Hfc) ■ z dS . (A.24) 

By analogy with the series expansion ( |A.18| ), an arbitrary function F(r^,2;) quadratically 
integrable on S 1 can be represented as a sum of the modal expansion *&(rt,z) in terms of 
eigenfunctions *&j-(rt) (the projection of F(rt,z) tangent to Hilbert space and convergent in 
mean) and the orthogonal complement C(rt,z): 

F(r t , z) = *(rt, z) + C(r t , z) = £ a fc (z) * fc (r t ) + C(r t , z) (A.25) 
where we have denoted 

F(r.,») = . = f^t i) . o(,„,) = 

\U(r t ,z)J \H a {r t ,z)J \H b {r t ,z)J 

and by analogy with Eq. ( |A.17| ) the orthogonal complement C satisfy the relation 

(* fc , C) = J ^l-t-CdS = J (E* k x H b + E b x H|.) ■ z dS = . (A.26) 

Eqs. ( |A25 ) and ( |A.26 ) allow the electromagnetic fields to be represented in the following 
form 

E(r t , z) = E a (r t , z) + E 6 (r t , z) = £ a k (z) E fc (r t ) + E 6 (r t , z) , (A.27) 

k 
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H(r t , z) = H a (r t , z) + H 6 (r t , z) = £ a k (z) H fc (r t ) + H 6 (r t) z) , (A.28) 

fe 

where the orthogonal complementary fields Ek(r£,z) and H&(rt, z) as well as the mode ampli- 
tudes a k (z) of the modal expansions 



B a (r t ,z) = J2a k (z)B k (r t ) and H a (r 4 , z) = ^ a k (z) ~H. k {r t ) 



(A.29) 



should be determined. The amplitudes afc(z), in principle, can be found from the equations 
similar to Eqs. ( |A.7| ) for the nonorthogonal basis or to Eq. ( A.15| ) for the orthogonal basis, with 
R k being given as follows 



&l ■ r • f ds 



J R fe = (* fc ,*) = (* fc ,F) = 

(EjxHlEx H* k ) • z dS 



in particular, by analogy with Eq. ( |A.15| ) 

R k _ (¥ fe ,F) 1 



(E^xHlExH^ z dS 



(A.30) 



(A.31) 



N k iVWs 

It is of great importance in electrodynamic applications that such a procedure of determining 
the mode amplitude a k (z) based on Eqs. (A. 7) or ( |A.15| ) allows us instead of the series expansion 
\l/(r t , z) in terms of eigenmodes to apply its finite sum of the nth order 



S n {v u z) = Y^a[ n \z)^ k (v t ) 



k=l 



like Eq. ( |A.2p , which yields the least mean-square error 



Si(r t ,z)-&(r t ,z) 



S n (r t ,z) - *(r t ,z) 



dS. 



(A.32) 



(A.33) 



analogously to Eq. ( |A,4j ). 

The above general reasoning concerning the convergence in mean, completeness, and orthog- 
onality properties of base functions can be extended to the electrodynamic basis of eigenfunc- 
tions so that, in particular, the generalized Parseval identity (All) and Bessel inequality ( A.19Q 
take the following form 



= ||*| 



*t r * ds 



is 

oo oo 



(E* a xH a + E a xH* a )-z Q dS=J2J2 N ki4 a i OT =E Ar fcl«fc! 2 ( A - 34 ) 

5 k=l 1=1 k=l 



and 



(F,F) = ||* + C| 



Ft • r • F dS 



„ oo oo oo 

/ (E* xH + Ex U*)-z dS> J2J2 N kia* k ai or > £ N k \ 

S k=l 1=1 k=l 



a k \ 



(A.35) 
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where the last single sums correspond to the orthogonal basis. 

As noted above, the orthogonality property of a basis is not mandatory but its existence 
facilitates the determination of the expansion coefficients at(z). Such a property is inherent in 
lossless physical systems, whereas losses destroy the "pure" orthogonality and convert it into 
the so-called quasi- orthogonality (see Sec. 3.1). 

B. DIRECT DERIVATION OF THE EQUATIONS OF MODE 
EXCITATION FROM MAXWELL'S EQUATIONS 



Starting point to derive the equation of mode excitation is formulas ( |4 . 1 2| ) and ( [4.13| ) which 
are a result of transforming Maxwell's equations Q4.4j ) and ( |4.5| ) inside the source region. Let 
us rewrite Eqs. (4.12) and ( 4.1 3| ) for transverse components: 



dAi 

^-p(z xH|) = -VxH 6 + iuPu + J 



V x Eb - iujfioM.bt — J 



bt > 



dz 



(B.l) 
(B.2) 



Here, in accordance with Eq. (2.2), the orthogonal complements for the polarization P b and 
magnetization M& are defined as 



P 6 = D 6 - e E b and /x M b = B b - /u H & 



(B.3) 



so that, as follows from Eqs. ( 4.14 ), (|4.15 ), ( 4.17 ), and ( 4.1S| ), their longitudinal components 
contribute to the complementary fields: 

1 



E b = zq E b 



Hfc = Zq H b 



iu>eo 
1 



;jf 2 + iuP bz ), 



-(J£ + iLJfi M, 



(B.4) 
(B.5) 



If we scalar-multiply Eqs. (|B.1| ) and (|B.2j ) by and — Et, respectively, and add the results, 
then after integrating over 5 we obtain 

dAi 

i 



J2^-J s (E* k x H z + E z x H* k ) -z dS = J s (E* k ■ V x H 6 - H£ • V x E 6 ) dS - 



{ (J e bt ■ E* kt + • H^) dS-iu f (P w • E^ + noMu ■ U* kt ) dS . (B.6) 
Js Js 



Now it is necessary to transform the first integral in the right-hand side of Eq. ( |B.6| ). The 
terms of its integrand can be rearranged as follows 

E* k ■ V x H 6 = z • (E* k x VH b ) = H b ■ V x E* k - z • V x (E* k H b ) = 
= iuno (H* k + M* k ) ■ H b - z • V x (El H b ) , 
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H£ • V x E b = z • (H£ x VE b ) = E b ■ V x H£ - z • V x (H* k E b ) = 
= -iu(e E* k + F* k )-E b - zq-Vx (H%E b ) , 

where in the last equalities for the fcth mode we have used Eqs. (|5.9| ) and ( |5.10|) . Then the first 
integral in Eq. ( B.6 ) turns into the sum of three integrals: 

/ (B* k -V xU b -Ul-V xE b )dS= f V x (H* k E b -E* k H b )-z dS + 

J S b J S b 

+ iu [ (E* k • e E b + U* k • Mo H 6 ) dS + icu [ (P* k • E 6 + ^M* k • H 6 ) dS . (B.7) 

J Sh J Sh 



The first integral in the right-hand side of Eq. ( [B.7D is transformed by using the Stokes 
theorem [42] into the following form 

/ V x (U* k E b -E* k H b )-z dS = I (H* k E b -E* k H b )-rdl = 

JS h JL b 

n b x H b ) ■ E£ — (n b x E b ) ■ H£ 



dl 



(Bi 



where n b and t = zq x n b are the unit vectors, respectively, normal (outward) and tangential 
to the contour L b bounding the bulk current area S b . 

The second integral in the right-hand side of Eq. flB.7| ) is rearranged by using Eqs. ( |B.4D 
and ( |B.5| ) to the following form 

iu [ (E* k • e E 6 + H£ • Mo H 6 ) dS = 

Js b 

= - f [J e bz E* kz + J% H* kz ) dS -iu f (P bz E* kz + ^M bz H* kz ) dS . (B.9) 



After inserting Eqs. ( |B.8[) and (|B.9| ) into Eq. (B/7) we obtain 

f (El • V x H 6 - H£ • V x E 6 ) dS=<f \(n b x H 6 ) • E* k - (n 6 x E 6 ) • H* k 

- [ (Jt z E* kz + J% H* kz ) dS-iuf (P bz E* kz + ^M bz H* kz ) dS + 
J s b J s b 

+ iu [ (P* k ■ E b + n M* k ■ H b ) dS . (B.10) 

Js b 

The first integral in the right-hand side of Eq. ( B.10| ) involves the effective surface currents 
n b x H b and J™ e/ = n fe x E b defined by Eqs. ( [4.271) and ( [4.281) . With allowing for this 
and employing the expression for the normalizing coefficient 



N kl 



[ (ElxHi + EixHD-zodS , 

Js 
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the substitution of Eq. ( |B.1C| ) into Eq. ( [B.6D yields 
dAi 

l 



N kl ~T~ = - I ( J ft ■ E fc + J™ • H D dS - I (J e s<ef • E* k + 3™ ef ■ H* k ) dl - 



ID 



(P 6 -K-K- E 6 ) + (wM 6 • H^, - Mo M£ • H 6 



(B.ll) 



The last integral in the right-hand side of Eq. ( B.ll ) vanishes because of 
(Pft • B* k - P* k • Eft) + ( M Mft • H* k - voMt • H ( 
(Dft • B% - D* k ■ Eft) + (Bft • U* k - B* k ■ Hft) 
(e-et) :E 6 E^ + (ft -ft t) : U b U* k + (|-C f ) :H b E^ + (C ~V) : E 6 H* fe 







where the constitutive relations Q2.8[ ), p.9| ), and ( p. 13 ) have been used for lossless bianisotropic 
media. 

For the most general case of the reactive kth. mode from the orthonormalization rela- 
tion ( p. 24 ) we have N k i = N k 5 kl . Then formula ( [B.ll ) finally gives the excitation equation for 
the kth mode: 

-L / (Jt-E* k + 3?-IS* k )dS- 

JS h 



dz 



N k J Sb 

^J Lh (KefK+ Refill) dl- 



(B.12) 



This formula is in agreement with the similar equation (5.12) obtained by another method, 
not counting the absence of the actual surface currents which can be considered as enclosed 
implicitly into the bulk currents. 
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